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Introduction 



"What lies beyond the Standard Model ?" is the question which intrigues many particle 
physicists today. It is surprising that such a question should arise when most of the existing 
data agree very well with the predictions of the Standard Model. However, there are strong 
theoretical arguments suggesting existence of new physics beyond the Standard Model. 
Independent of these arguments, there are some recent experimental results which may be 
considered as the first indications of some new physics. Put together, these arguments 
constitute a 'strong evidence' for the presence of physics beyond the Standard Model and 
thus the above question. 

At present, the question has not been answered. Future experiments like LHC (Large 
Hadron Collider) are likely to provide an answer. To satisfy some of the theoretical ar- 
guments against it, the Standard Model is usually extended by assuming an additional 
symmetry called super symmetry. One of the most popular of such models is the Minimal 
Supersymmetric Standard Model (MSSM). The popularity of this model lies in the fact that 
it not only satisfies one of the major theoretical prejudices, namely the hierarchy problem, 
but is also testable in the future colliders. On the other hand, the experimental 'evidence' 
for physics beyond the Standard Model comes from the recent results of the atmospheric 
neutrino experiments which indicate that the neutrinos may have small masses and non-zero 
mixing among them. The Standard Model which has no provision for neutrino masses thus 
has to be extended in some sector viz., with additional symmetries, additional particles etc. 
These extensions are achieved in many cases, independent of the theoretical prejudices, thus 
leading to a different set of extensions of the Standard Model. 

In this thesis, we assume supersymmetry as the physics beyond the standard model. We 
then study the generation of neutrino masses within supersymmetric standard models and 
implications from the results of the solar and atmospheric neutrino experiments on these 
models. But, before proceeding further, we try to motivate our work in this chapter, along 
with a summary of salient features of the MSSM. 

1.1 The Standard Model 

As a starting point, we here briefly review the salient features of the Standard Model. 
The Standard Model (SM) Q is a spontaneously broken Yang-Mills quantum field theory 
describing the strong and electroweak interactions. The theoretical assumption on which 
the Standard Model rests on is the principle of local gauge invariance with the gauge group 
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given by, 

GsM = SU{?>)c X SU{2)l X U{1)y (1.1) 

The particle spectrum and their transformation properties under these gauge groups are 
given as, 

Qi = ( "^^M ~ fs, 2, l\ Ui = URi ~ (3, 1, I 



A = c^fli ~ (3, 1, - ^ 
= I ''"'O ~ fl' 2, - I) Ei = CR, ~ (1, 1, - 1) 



2 



(1.2) 



In the above i = 1,2,3 stands for the generation index. Qi represents the left handed quark 
doublets, Li represents left handed lepton doublet, Ui, Di, Ei represent right handed up- 
quark, down-quark and charged lepton singlets respectively. The numbers in the parenthesis 
represent the transformation properties of the particles under Gsm in the order given in 
eq.([l.l|). In addition to the fermion spectra represented above, there is also a fundamental 
scalar called Higgs whose transformation properties are given as, 

-(^;)^(..i) 

There are also gauge boson fields which enter the spectrum through the requirement of 
local gauge invariance. The total lagrangian with the above particle spectrum and gauge 
group can be represented as, 

C-SM = + CyM + Cs + Cyuk- (1-4) 
The fermion part, £i? is given as, 

Cf = i^-i^V^^ (1.5) 

with 

^ = {Qi Ui, Di, Li, Ei) (1.6) 
where T>^ represents the covariant derivative of the field given as, 

V^ = dii- igsG^X^ - i^Wy - ig'B^Y (1.7) 

Here A = 1,..,8 with representing the SU{2>)c gauge bosons, / = 1,2,3 with 
representing the SU{2)l gauge bosons. The U{1)y gauge field is represented by B^. The 
self interactions of the gauge fields are given by, 

Cym = -\g^''^G% - -^^""'Wl, - -B^^'B,, (1.8) 
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with 

G'^u = dfj^G^ - dyG^ + Qs JabcG^G'^ 

B^u = d^B,-d,B^, (1.9) 

where fABC(iJK) represent the structure constants of the SU{3){SU{2)) group. 
The scalar part of the lagrangian is given by, 

Cs = (V^H)^ V^H - V{H), (1.10) 

where 

V{H) = n^H^H + A {h^H^ ^ (1-11) 
and finally the Yukawa part is given by, 

J^yuk = KjQiUjH + hfjQiDjH + hfjLiEjH + H.c (1.12) 

where H = ia^H*. The symmetry Gsm is not realised in nature. It has to be broken to 
GsM — > SU (3)c X U{l)em- This is achieved in the Standard Model in a renormalisable way 
by spontaneous symmetry breaking. The Higgs field attains a vacuum expectation value 
{vev) in this mechanism called the 'Higgs mechanism'. Fermions (except neutrinos) and 
gauge bosons attain mass through this mechanism. Only one degree of freedom remains for 
the Higgs field, the Higgs boson. 

The Standard Model is renormalisable and anomaly free |Q] . The fermion spectra 
of the SM has been completely discovered with the recent discovery of Tau-neutrino 
(i't)- In addition to verifying the predicted neutral current interactions, experiments are 
now able to probe higher order radiative corrections with precision measurements The 
results from these experiments agree with the predictions of the SM very well There 
has been no indication so far for deviations from the Standard Model physics Q in this case. 
On the other hand, the scalar sector of the SM is not very well known. The Higgs boson is 
yet to be discovered. Precision observables suggest a light Higgs scalar ~ 60 — 200 GeV , 
whereas the present experimental bounds put it within 106 GeV < uih < 235 GeV Q. 
The first number comes from the non observation of e~^e~ ZH at LEPII, whereas the 
upper bound comes from a global fit to the data [^. Thus we see that the Standard Model 
stands on a firm footing both theoretically and experimentally, except for the scalar sector 
which is yet to be discovered. Incidentally, it is the existence of a fundamental scalar in the 
Standard Model that gives rise to most of the conceptual problems with the SM, when SM 
is being incorporated in a bigger theory. But before that we will see the need for extending 
the Standard Model. 

In spite of the success of the Standard Model in explaining most of the observed phe- 
nomena, there are several features of the SM which still need to be understood. The SM 
contains about 20 parameters which are fixed by the experiments. Thirteen of them are 

^Please see |7|, for a recent fit to precision electroweak observables and their comparison to SM predictions. 
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related to the masses of the quarks and leptons ( nine fermion masses, three mixing angles 
and one phase of the CKM matrix). All these masses which are generated by the Yukawa 
part of the lagrangian, Cyuk (eq. |1.12| ) are hierarchical in nature in the generation space. 
At present, we have no understanding of the origin of this hierarchical nature. Similarly, we 
have no understanding of why there are only three generations within the SM. All these are 
collectively known as the 'flavour problem' of the SM. As a solution to this problem, one 
typically extends the SM with an additional 'flavour symmetry'. The quantum numbers of 
the SM fields under this symmetry would now decide the magnitude of the masses these 
particles take. 

Whereas the above arguments require an extension of the SM for a deeper understanding 
of the Standard Model, recent results from the atmospheric neutrino experiments make 
this requirement almost a necessity. The atmospheric neutrino experiments measure the 
flux of the electron and muon neutrinos produced in the atmosphere due to cosmic ray 
collisions with the air molecules. These flux measurements (or rather the ratio of the 
fluxes ) are anomalous with respect to the predicted flux from various calculations. This 
reduction of the measured flux is generally understood in terms of 'neutrino oscillations', 
which require neutrinos to have small masses. Recently, super-Kamiokande experiment 
has reported evidence for the existence of neutrino oscillations |1C]. If taken seriously, 
these results would imply existence of neutrino masses. In addition to these signals from 
the atmospheric neutrinos, there are other neutrino experiments like the solar neutrino 
experiments which also suggest neutrinos to have small masses. As we have seen earlier, the 
SM does not have right handed neutrinos in its spectrum eq. (|1.2| ), thus denying neutrinos 
any mass through Cy^k Q- Thus one has to extend the SM in some sector (symmetries, 
particles or both) to generate mass for the neutrinos. One of the most standard methods to 
generate neutrino mass is to add right-handed neutrinos in to the Standard Model particle 
spectrum. Gauge invariance would allow not only Dirac mass terms but also Majorana mass 
terms for the right handed neutrino fields. The interplay between these terms gives rise to 
a 'see-saw' mechanism, in which the left handed neutrinos attain small majorana masses. 

Thus we have seen that the SM has to be extended in order to satisfy requirements 
coming from neutrino experiments as well as for a deeper theoretical understanding. There 
is also a search for a Grand Unified Theory (GUT) |11] which would unify the strong and 
the electroweak forces with a single coupling constant. It is interesting to note that in- 
addition to unification of forces, neutrino masses can also be naturally incorporated within 
the Grand Unified theories |]l^] through the 'see-saw' mechanism discussed above. The 
Standard Model would then be an effective theory valid up to some scale relevant to the 
weak scale. The scale at which the GUT theories would take over from the SM is estimated 
from the life time of proton and 'Renormalisation Group' running of the three coupling 
constants to be ~ 10^^ GeV [|l3| |^. Inclusion of the gravitational forces would push this 
scale further up to the Planck scale, Mpianck ~ 10^^ GeV. This would imply that the SM 
has to be valid right from the weak scale up to the GUT-scale typically, fourteen orders 
of magnitude. This would create a conceptual problem within these theories called the 

^Neutrinos cannot attain majorana masses either, as lepton number is 'accidentally' conserved in the SM. 
^ eV-scale neutrino masses in some GUT theories would also require a scale of same order. 
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hierarchy problem. 

The hierarchy problem arises due to the existence of a fundamental scalar within the 
Standard Model A typical Grand Unified Theory uses a single gauge group and a 

single coupling constant to describe the physics at the GUT scale, Mqut ~ 10^^ GeV. This 
symmetry is then spontaneously broken at those scales to the symmetry of the Standard 
Model, GsM given in eq. (|l.l|) . As a consequence of this breaking, the gauge bosons medi- 
ating the grand unified interactions acquire heavy masses of the order of Mqut and thus 
would not be relevant for the physics at the weak scale. The gauge bosons of the Standard 
Model remain massless as Gsm remains unbroken. The fermion masses also would not be 
affected by this breaking as they are protected by chiral symmetries. But, the scalar mass is 
unprotected by any symmetry. There is no reason to assume the scalar would not gain mass 
of the order of Mqut due to the spontaneous breaking of Gqut — > Gsm 0|- Such a large 
mass for the scalar is disastrous phenomenologically as the observed masses of the gauge 
bosons typically require a scalar mass of O(IOO) GeV. To solve this 'hierarchy' problem, one 
can either bring down the scale of the unification or introduce an extra symmetry to protect 
the Higgs mass. The former would require introduction of extra space time dimensions |l5[ 
whereas the symmetries which would protect the Higgs mass are called super symmetries. 
There is also one more approach in the literature which assumes that the Higgs boson is a 
composite particle of fermions thus ruling out the existence of a fundamental scalar in the 
theory ||l6|| . 

In the present thesis, we follow the approach of supersymmetry. Supersymmetry in ad- 
dition to protecting the Higgs mass from attaining large values has other attractive features. 
A minimal super symmetric version of the standard model can be built which is predictable 
and testable at present and future colliders. The unification of gauge coupling constants at 
high scales is exact in this case and so. Grand Unified Theories also prefer the existence of 
supersymmetry. More fundamental theories like string theories may also require supersym- 



metry to be present as low energy effective theories [17|. Thus super symmetric standard 
models make an attractive framework as the physics beyond the Standard Model. Below, 
we review some salient features of the supersymmetric version of the Standard Model. 

1.2 The Minimal Supersymmetric Standard Model 

Supersymmetry is a symmetry that transforms a fermion in to a boson and vice versa. 
To understand how it protects the Higgs mass, let us consider the hierarchy problem once 
again. The Higgs mass enters as a bare mass parameter in the Standard Model lagrangian, 
eq.( |l.lO ). Contributions from the self energy diagrams of the Higgs are quadratically di- 



vergent pushing the Higgs mass up to cut-off scale |18|. In the absence of any other new 
physics at intermediate energies, the cut-off scale is typically Mgut or Mpianck- Cancella- 
tion of these divergences with the bare mass parameter would require fine-tuning of order 
one part in 10~^^ rendering the theory 'unnatural' [19|. On the other hand, if one has an 



additional contribution from a fermionic loop, with the fermion carrying the same mass as 
the scalar, the contribution from this diagram would now cancel the quadratically diver- 
gent part, with the total contribution now being only logarithmically divergent [KG]. This 
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would require a symmetry which would relate a fermion and a boson with same mass. Such 
symmetries are known as supersymmetries. 

Supersymmetries were first introduced in the context of string theories by Ramond 



l|. In quantum field theories, this symmetry is realised through fermionic generators [22| 



thus escaping the no-go theorems of Coleman and Mandula |23]. The simplest Lagrangian 



realising this symmetry in four dimensions was built by Wess and Zumino ||2j| which contains 
a spin ^ fermion and a scalar. To build interaction lagrangians one normally resorts to 



the superfield formalism of Salam and Strathdee |25], as superfields simplify addition and 
multiplication of the representations ||2^. It should be noted however that the component 
fields may always be recovered from superfields by a power series expansion in grassman 
variable. 

A chiral superfield contains a weyl fermion, a scalar and and an auxiliary scalar field 
generally denoted by F. A vector superfield contains a spin 1 boson, a spin 1/2 fermion 
and an auxiliary scalar field called D. A minimal supersymmetric extension of the Standard 
Model 1 27 1 is built by replacing every standard model matter field by a chiral superfield 
and every vector field by a vector superfield. Thus the existing particle spectrum of the 
Standard Model is doubled. The particle spectrum of the MSSM and their transformation 
properties under Gsm is given by, 



(1.13) 

The scalar partners of the quarks and the leptons are typically named as 's'quarks and 
's'leptons. For example, the scalar partner of the top quark is known as the 'stop'. In the 
above, these are represented by a 'tilde' on their standard counterparts. A distinct feature 
of the supersymmetric standard models is that they require two Higgs fields with opposite 
hypercharges to cancel the anomalies and to give mass to both up-type and the down-type 
quarks ^ These transform as 




*The same notation is followed in the entire thesis. 

^The hermitian conjugate of the superfield is not allowed into the superpotential due to its holomorphic 



nature. 
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The fermionic parts of the Higgs superfields are known as Higgsmos. The gauge bosons and 
their fermionic partners, gauginos are represented as, 



: ( G^'^ 
: ( 

Vy : [b^^ b) (1.15) 

In the above, the fermionic (susy) counter parts of the Higgs and the gauge bosons have 
been represented by a 'tilde' on their bosonic notations. The total Lagrangian of the MSSM 
is of the form : 

Cmssm = j d6^ W{^) + j de^ dO^ $1 e^^ + j dO^ W^W^. (1.16) 

In the above, VF(<I>) is a function of the chiral superfields called the superpotential, $ 
representing a generic superfield. The renormalisable superpotential is of dimension three 
or less. For the MSSM, the superpotential invariant under Gsm is given as, 



W = Wi + W2, (1.17) 



where 



Wi = h'^jQiU^H2 + hfjQiD^Hi + htjLiE^Hi + fiHiH2 (1.18) 
W2 = e,L,H2 + \ijkULjEl + \[^kUQjDl + >;ijkUiD''pl. (1.19) 

The second term in the RHS of eq.( |l.l6| ) represents the matter-gauge boson coupling, 
with V = {yf^-,V^,V^^ with the appropriate coupling constants. The last term in the 
RHS of the eq.( |1.16| ) represents the self interaction of the gauge fields. The Lagrangian 
in component form can be found by expanding the superfields as noted above. The scalar 
potential can be derived by eliminating the auxiliary fields D and F which appear in the 
definitions of the vector superfield and chiral superfield respectively. 

1.3 R-parity 

In the above, we have seen the minimal supersymmetric version of the Standard Model. 
Comparing eq. (|1.17 ) with the Standard Model Lagrangian, we see that in additional to 



doubling of the particle spectrum, new interactions which violate baryon number and in- 
dividual lepton number are now allowed in the lagrangian. This is because the Higgs 
superfield. Hi and the lepton superfields Lj have the same quantum numbers under Gsm- 
Though they carry the same quantum numbers in the SM also, these interactions are not 
possible as the fermions and the Higgs transform as different representations under Lorentz 
Group. Within the MSSM the distinction is lost since the superfields corresponding to the 
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leptons and the Higgs transform identically under supersymmetry and gauge symmetry and 
these interactions appear as they are gauge invariant. 

The first three terms of the second part of the superpotential W2 (eq.( |l.l^ ), are lepton 
number violating whereas the last term is baryon number violating. The simultaneous 
presence of both these interactions can lead to proton decay for example, through a squark 
exchange. Stringent limits can be placed on the products of these couplings from the life 
time of proton To avoid proton decay, one can either remove both these couplings or 
assume one set (either baryon number violating or lepton number violating couplings) to be 
zero. The former is normally arrived at by imposing a discrete symmetry on the lagrangian 
called R-parity. R-parity has been originally introduced as a discrete R-symmetry ^ by 



Ferrar and Fayet |29] and then later realised to be of the following form by Ferrar and 
Weinberg pO| ]: 

where B and L represent the Baryon and Lepton number respectively and F represents the 
Fermion parity given as -1 for fermions and + 1 for bosons. Under R-parity the transfor- 
mation properties of various superfields can be summarised as: 



9 -r 

{Qi,Uf,DlLi,Ef} ^ -{Qi,U^,D^,U,En 

{Hi,H2} ^ {Hi,H2} (1.21) 

Imposing R-parity has an advantage that it provides a natural candidate for dark matter. 
This can be seen by observing that R-parity distinguishes a particle from its superpartner. 



Thus the lightest supersymmetric particle (LSP) cannot decay and remains stable |31]. R- 
parity has also been motivated as a remnant symmetry in a large class of supersymmetric 
theories, especially in theories with an additional U{1) symmetry. It has also been pointed 
out that though R-parity constraints baryon and lepton number violating couplings of 
dimension four, dim 6 operators which violate baryon and lepton numbers are still allowed 



by R-parity [32| 



1.4 Supersymmetry breaking 

The MSSM lagrangian built so far is invariant under supersymmetry. Supersymmetry 
breaking has to be incorporated in the MSSM to make it realistic. In a general lagrangian, 
supersymmetry can be broken spontaneously if the auxiliary fields D or F appearing in 
the definitions of the chiral and vector superfields respectively attain a vacuum expectation 
value (vev). Incorporation of this kind of breaking in MSSM using the MSSM superfields 
leads to phenomenologically inconsistent scenarios, like for example existence of a very 

®R-symmetries are symmetries under which the 9 parameter transform non-trivially. 
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light scalar Q An alternative approach to incorporate supersymmetry breaking is using the 
'hidden-sector' scenarios. 

The hidden-sector is a sector of superfields which do not carry any quantum numbers 
under ordinary gauge interactions, i.e, Gsm- The only interactions they share with the 
visible (MSSM) sector superfields is through gravity. Supersymmetry can then be broken 
spontaneously in the hidden-sector and this information is then passed on to the MSSM 
sector through gravitational interactions. Since gravitational interactions play an important 
role only at very high energies, Mp ~ 0(10^^) GeV, the breaking information is passed on 
to the visible sector only at those scales. The end effect of this mechanism is that explicitly 
supersymmetry breaking terms are now added in to the lagrangian. These 'soft' terms 
do not reintroduce quadratic divergences back into the theory. The strength of the soft 
terms is characterised by, Mg/Mpianck^ where Ms is the typical scale of supersymmetry 



breaking. These masses can be comparable to weak scale for Ms ~ 10 GeV [34]. The 'soft' 



super symmetric breaking part of the lagrangian can be classified to contain the following 



terms Pq] : 



• a) Mass terms for the gauginos, Mi, M2, M^. 

• b) Mass terms for the scalar particles, m^, with (pi representing the scalar partner of 
chiral superfields of the MSSM. 

• c) Trilinear scalar interactions, Aijk(j)i(j)j(f)j. corresponding to the cubic terms in the 
superpotential. 

• d) Bilinear scalar interactions, Bij(l)i(j)j corresponding to the bilinear terms in the 
superpotential. 

The total MSSM lagrangian is thus equal to 

C-total = C,MSSM + C,soft (1-22) 

with Cmssm given in eg. ([1.161). 



1.4.1 Universality and CMSSM 

The above mechanism of supersymmetry breaking is called minimal supergravity (mSUGRA) 
inspired supersymmetry breaking. As we have seen above this type of breaking introduces 
several new soft parameters in to the theory. Typically the number of these parameters is 
large ~ 105. Moreover, large flavour changing neutral currents (FCNC) are also introduced 
by this kind of breaking |3(:]. To reduce the number of parameters as well as remove the 



large FCNC contributing terms, an ansatz is usually followed at the high scale where the 
soft terms are introduced in to the theory. This ansatz is called universality and it assumes 
the following: 



''This can be seen from the mass sum rules which appear as a consequence of spontaneous supersymmetry 
breaking |M. 
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• All the gaugino mass terms are equal at the high scale. 

Ml = M2 = M3 = M 

• All the scalar mass terms at the high scale are equal. 

2 2 

• All the trilinear scalar interactions are equal at the high scale. 

• All bilinear scalar interactions are equal at the high scale. 



It is possible to construct supergravity models which can give rise to such kind of strong 
universality [^]. This approximation now drastically reduces the number of parameters of 
the theory to about five, mQ,M (or equivalently M2), ratio of the vevs of the two Higgs, 
tan/3, A, B. Thus, these models are also known as 'Constrained' MSSM [^] in literature. 
The low energy mass spectrum of the soft terms is now determined by the Renormalisation 
Group scaling of those parameters. The supersymmetric mass spectrum of these models has 
been extensively studied in literature [^]. The Lightest Supersymmetric Particle (LSP) is 
mostly a neutralino in this case. 

1.4.2 Gauge Mediated Supersymmetry breaking 

In the above we have seen that supersymmetry is broken at a high scale and is communicated 
through gravity to the normal particle sector. It induces large FCNC's and a large number 
of parameters which can be corrected by assuming a 'strong' universality at the Mqut- 
An alternative mechanism has been proposed which tries to avoid these problems in a 
natural way. The key idea is to use gauge interactions instead of gravity to mediate the 
supersymmetry breaking from the hidden (also called secluded sector sometimes) to the 
visible MSSM sector [^^. In this case supersymmetry breaking can be communicated at 
much lower energies ~ 100 TeV. 

A typical model would contain a susy breaking sector called 'messenger sector' which 
contains a set of superfields transforming under a gauge group which 'contains' Gsm- Su- 
persymmetry is broken spontaneously in this sector and this breaking information is passed 
on to the ordinary sector through gauge bosons and their fermionic partners in loops. The 
end-effect of this mechanism also is to add the soft terms in to the lagrangian. But now 
these soft terms are flavour diagonal as they are generated by gauge interactions. The soft 
terms at the messenger scale also have simple expressions in terms of the susy breaking 
parameters. In addition, in some models of gauge mediated supersymmetry breaking, only 
one parameter can essentially determine the entire soft spectrum. 
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In a similar manner as in the above, the low energy susy spectrum is determined by 
the RG scaling of the soft parameters. But now the high scale is around 100 TeV instead 
of Mqut as in the previous case. The mass spectrum of these models has been studied in 
many papers |41]. The lightest supersymmetric particle in this case is mostly the gravitino 
in contrast to the mSUGRA case. 



1.4.3 SU{2) X [/(I) breaking 

As we have seen earlier, the supersymmetric version of the Standard Model is a two Higgs 
doublet model. A consistent incorporation of the SU{2)l x U{1)y breaking puts constraints 
relating various parameters of the model. To see this, consider the neutral Higgs part of 
the scalar potential. It is given as 



Vscalar = {mj,^ + fi^)\Hf\^ + {mj,^ + fM^)\H^\^ - {Bf,^iH^ + H .c) 

+ lia' + 9"){Hf - Hff + . . . , (1.23) 

where Hi,H2 stand for the neutral Higgs scalars and we have parameterised the bilinear 
soft term B = B^fj,. The existence of a minima for the above potential would require 
at least one of the Higgs mass squared to be negative. In both gravity mediated as well 
as gauge mediated supersymmetry breaking models, such a condition at low energies can 
be naturally incorporated. The soft parameters which appear in the above potential are 
generally positive at the high scale. But radiative corrections significantly modify the low 
scale values of these parameters, making one of the Higgs mass to be negative at the weak 
scale. This mechanism is called radiative electroweak symmetry breaking. In addition to 
ensuring the existence of a minima, one would also require that the minima should be able 
to reproduce the standard model relations i.e, correct gauge boson masses. This would 
give rise to constraints on the parameters known as the minimisation conditions. These are 
given as 



2^ - t^ 



o^i = \.^2 a \ ' -f^' (1-24) 



Sin2/? = , (1.25) 

where tan (3 = V2/V1 is the ratio of the vacuum expectation values of H2 and respec- 
tively 0. In addition to the above conditions care should also be taken such that charge and 
color breaking minima are absent. 

^The above minimisation conditions are given for the 'tree level' potential only. The minimisation con- 
ditions for the one-loop effective potential are given in |42l. 
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1.5 Motivation and Outline of thesis 

So far in this chapter we have seen the reasons for beheving in the existence of physics 
beyond the Standard Model. A strong signal comes from the recent neutrino oscillation 
experiments. On the other hand, extension scenarios with unification ideas are generally 
plagued by the hierarchy problem. A solution to this problem is assuming supersymmetry 
just above the weak scale. Supersymmetric standard models are built which can incorporate 
SU (2) X U{1) breaking naturally. The question then remains is whether one can incorporate 
neutrino masses and mixing within these theories in a natural way. 

Supersymmetric Standard Model, unlike the Standard electroweak model has a natural 
source of lepton number violation. Since, there is also baryon number violation in these 
theories, one imposes R-parity to remove both these set of couplings. But R-parity is 
ad-hoc. One can always assume symmetries other than R-parity like for example baryon 
parity ||4^ which can remove the baryon number violating couplings, leaving us with lepton 
number violating couplings only. In the presence of these lepton number violating couplings, 
neutrinos attain major ana masses f\. The generation of the neutrino masses in these theories 
and whether they are of the correct order to satisfy the solar and atmospheric neutrino 
masses is the main subject of this thesis. In addition to neutrino masses R-parity violating 
theories have a different set of experimental signatures in complete contrast to the standard 



MSSM signatures ||4J. We also study some such experimental signatures in the context of 
some recent experimental results of HERA detector. 

The outline of the thesis is as follows. In chapter 2, we discuss the various neutrino 
experiments and some standard neutrino mass models. In chapter 3, we study the standard 
Renormalisation group analysis of the MSSM and then try to understand the effect of 
RG analysis on the neutrino mass spectrum. In chapter 4, we discuss a specific model of 
R-violation namely bilinear R-violation within the framework of gauge mediated models 
supersymmetry breaking. In chapter 5, we discuss models with trilinear lepton number 
violation and the feasibility of simultaneous solutions to solar and atmospheric neutrino 
problems. In chapter 6, we derive bounds on the trilinear A'^^^^ couplings from neutrino 
mass and study its implications for the HERA anomalies. We end with conclusions and 
future outlook in chapter 7. 



®Even with R-conservation, higher dimensional R-parity violating operators would give rise to small 
neutrino masses fB2|. 



References 



[1] The original papers of the electroweak theory are given by S. L. Glashow Nucl. Phys. 
22 (1961) 579; S. Weinberg Phys. Rev. Lett. 19 (1967) 1264 ; A. Salam, in Elementary 
Particle Physics ed. N. Svarthom pg. 367; The strong interactions were first described 
by, D. Gross and F. Wilczek Phys. Rev. Lett. 30 (1973) 1343; D. Pohtzer, Phys. Rev. 
Lett. 30 (1973) 1346. The Standard Model has also been reviewed in S. Abers and B. W. 
Lee, Phys. Rep. 9 (1973) 1; E. D. Commins and P. H. Bucksbaum, Weak Interactions 
of Leptons and Quarks, Cambridge University Press (1983) Chapter 2. Recent reviews 
include, W. Herrero, Lectures presented at NATO-ASI school, hep-ph/9812242. 

[2] G. 't Hooft, Nucl. Phys. B35 (1971) 167. 

[3] R. Jackiw in Lectures in Current Algebra and its Applications, Princeton University 
Press, Princeton (1972); Lewis H. Ryder Quantum Field Theory Second Edition, 
Cambridge University Press, (1996), Chap. 9.10, pg. 366; S. Weinberg, The Quantum 
Theory of Fields, Cambridge University Press (1996), Chapter 22. 

[4] D. E. Groom et al, The Review of Particle Physics, The European Physical Journal 
C15, 1 (2000); The Particle Data Group, http://pdg.lbl.gov/. 

[5] P. Langacker, Lectures presented at Theoretical Advanced Training Institute in Ele- 
mentary Particle Physics, (TASI-92), Boulder, U. S. A, hcp-ph/9303304 ; W. Hollik, 
Talk given at 5th Hellenic School and Workshops on Elementary Particle Physics, 
(1995), Corfu, Greece, hep-ph/9602380 ; See also, G. Azuelos et al hep-ph/0003275 
for prospects at LHC. 

[6] The status of the Standard Model has been reviewed by many authors recently. Here we 
state some of them. W. Hollik, talk given at ICHEP '98, hep-ph/9811313; J. Erler and 
P. Langacker, talk presented at WEIN '98, hep-ph/9809352 ; J. L. Hewett, Lectures 
given at TASI-97, hep-ph/9810316. 

[7] F. Zwirner in , Particle Physics : Ideas and Recent Developments, Ed. J. J. Aubert, R. 
Gastmans and J. M. Gerard, Kluwer Academic Publishers, 2000. (Lectures presented 
at Cargese Summer Institute, Cargese, Prance, 1999.) 

[8] For a recent review, please see, W. J. Marciano, in the proceedings of 5th International 
Conference on Physics Potential and Development of Muon Colliders (MUMU 99), San 



13 



References 



14 



Francisco, California, (1999), hep-ph/0003181. Also see, J. Erler, hep-pli/0010153 and 
R. Barbieri and A. Strumia, liep-ph/0007265. 

[9] J. Erler, Talk presented at 17th International Workshop on Weak Interactions and 
Neutrinos (WIN 99), Cape Town, South Africa, Jan 1999, hep-ph/9904235. 

[10] Y. Fakuda et. al, super-Kamiokande Coll. Phys. Rev. Lett. 81 (1998) 1562. 

[11] R. N. Mohapatra, Supersymmetry and Grand Unification; Lectures given at ICTP 
Summer School in Particle Physics, Trieste, Italy, hep-ph/9911272 ; Lectures at TASI 
97, hcp-ph/9801235. For an introduction, please see G. G. Ross, Grand Unified Theo- 
ries, Benjamin/cummings ( 1984) ( Frontiers In Physics, 60). 

[12] M. Gell-man, P. Ramond and R. Slansky, in Supergravity, North Holland, Amsterdam, 
(1979). R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

[13] P. Langacker, Phys. Rep. 72 (1981) 185. 

[14] E. Gildener and S. Weinberg, Phys. Rev. D13 (1976) 3333. 

[15] I. Antoniadis, Phys. Lett. B246 (1990) 377; N. Arkani-Hamed, S. Dimopoulos and G. 
Dvah, Phys. Lett. B429 (1998) 263; I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos 
and G. Dvah, Phys. Lett. B436 (1998) 257. 

[16] K. Lane, Lectures given at LNF Spring School in Nuclear Subnuclear and Astroparticle 
Physics, Rome, Italy (2000), hep-ph/0007304; Lectures given at TASI-'93, Boulder, 
CO, hep-ph/9401324. Also see, K. Lane, hep-ph/9610463. 

[17] E. Kiritsis, Particle Physics : Ideas and Recent Developments, Ed. J. J. Aubert, R. 
Gastmans and J. M. Gerard, Kluwer Academic Publishers, 2000. (Lectures presented 
at Cargese Summer Institute, Cargese, France, 1999.); M. Dine, hep-ph/001035. 

[18] L. Susskind, Phys. Rev. D20 (1979) 2019 

[19] G. 't Hooft, Recent Developments in Gauge Theories, (1980) Plenum Press. Also see, K. 
Agashe, Ph. D. Thesis, hep-ph/9908491 chapter 1 for a nice discussion on naturalness 
in supersymmetric standard models. 

[20] Michael Dine, Lectures given at TASI-'98, hep-ph/9612389. 

[21] P. Ramond, Phys. Rev. D3 (1971) 2415; Quantum Field Theory: A modem Primer, 
Addison- Wesley (1989) (frontiers in physics, 74). 

[22] R. Haag, J. T. Lopuszanksi and M. Sohnius, Nucl. Phys. B88 (1975) 257; P. Fayet and 
S. Ferrara, Phys. Rep. 32 (1977) 249; S. Weinberg, The Quantum Theory of Fields 
Vol. Ill, Cambridge University Press (2000), chapter 1. 



References 



15 



[23] S. Coleman and J. Mandula, Phys. Rev. D3 (1967) 1251. See Also, S. Weinberg, 
The Quantum Theory of Fields Vol. Ill, Cambridge University Press, 2000, chapter 1, 
Appendix B. 



[24 
[25 
[26 

[27 

[28 
[29 
[30 
[31 

[32' 
[33 

[34 
[35; 
[36 

[37; 
[38 



J. Wess and B. Zumino, Nucl. Phys. B70 (1974) 39. 

A. Salam and J. Strathdee, Nucl. Phys. B76 (1974) 477. 

J. Wess and J. Bagger Supersymmetry and Supergravity, Princeton University Press, 
1983; H. P. Miles, Phys. Rep. 110 (1984) 1; J. Bagger, TASI-96, liep-ph/9604232 ; S. 
Weinberg Vol III, Chapter 26. 

H. E. Haber and G. L. Kane, Phys. Rep. 117 (1985) 75; C. Csaki, Mod. Phys. Lett. 
All (1996) 599, hep-ph/9606414 ; S. P. Martin, hep-ph/9709356. 

A. Yu. Smirnov and F. Vissani, Nucl. Phys. B460 (1996) 37. 

G. Ferrar and P. Fayet, Phys. Lett. 76B (1978) 575. 

G. Ferrar and S. Weinberg, Phys. Rev. D27 (1983) 1731. 

J. S. Hagehn, D. V. Nanopulos, K. A. Olive and M. Sredinicki, Nucl. Phys. B238 
(1984) 453; J. Ellis, T. Falk, K. A. Olive and M. Schmitt, Phys. Lett. B413 (1997) 
355. Also see, J. Ellis, astro-ph/9812211. 

S. Weinberg, Phys. Rev. D26 (1982) 287. 

S. Ferrara, L. Girardello and P. Fayet, Phys. Rev. D20 (1979) 403; P. Fayet, Phys. 
Lett. B84 (1979) 416. Also see, S. Dimopuolos and H. Georgi, Nucl. Phys. B193 
(1981) 150. 



H. P. Nilles in ref.|2|]. 

L. Girardello and M. T. Grisaru, Nucl. Phys. B194 (1982) 65. 

A. Masiero and L. Silvestrini, Lectures given at International School of Subnuclear 
Physics, 35th Course: Highlights: 50 Years Later, Erice, Italy, hep-ph/9711401. Also 
see, F. Gabbiani et al , Nucl. Phys. B477 (1996) 321. 

H. Murayama, Lectures presented at ICTP summer school 1999, hep-ph/0002232. 

A. Djouadi et. al, The minimal supersymmetric standard model: group summary re- 
port, hep-ph/9901246. 



[39] For reviews please see, M. Drees, Pramana, 45 (1995) S85, hep-ph/9407254; M. Drees 
and S. P. Martin, hep-ph/9504324; M. Drees, hep-ph/9611409 and references there in. 
See also for collider physics, S. Abdullin et. al, The susy working group: summary 
report, hep-ph/0005142. 



References 



16 



[40] M. Dine and A. E. Nelson, Phys. Rev. D48 (1993) 1277; M. Dine, A. Nelson and Y. 
Shirman, Phys. Rev. D51 (1995) 1362 ; M. Dine, A. Nelson, Y. Nir and Y. Shirman, 
Phys. Rev. D53 (1996) 2658. 

[41] R. Rattazzi and U. Sarid, Nucl. Phys. B501 (1997) 297; S. Dimopuolos, S. Thomas and 
J. D. Wells, Nucl. Phys. B488 (1997) 39. For a review please see, S. L. Dubovsky and 
D. S. Gorbunov and S. V. Troitsky, hep-ph/9905466; G. F. Giudice and R. Rattazzi, 
liep-pli/9801271. 

[42] C. L. Mouel and G. Moultaka, Talk given at Meeting on Tests of the Electroweak 
Symmetry Breaking, Ouranoupolis, Greece, May 1997, hep-ph/9708368; Also see, C. 
L. Mouel and G. Moultaka, Nucl. Phys. B518 (1998) 3. 

[43] L. E. Ibanez and G. G. Ross, Nucl. Phys. B368 (1992) 3. 

[44] H. Dreiner, hep-ph/9707435; G. Bhattacharyya, hep-ph/9709395; R. Barbier et. al, 
Summary of the working group on the R-parity violation, hep-ph/9810232. 



Chapter 2 



Neutrino Anomalies and Mass Models 



Neutrinos are one of the most abundant Standard Model particles in the universe. In 
addition to the various natural radioactive sources, neutrinos are produced in the hydrogen 
burning process in the stars as well as when a star dies in supernova explosions. They are 
also produced when energetic cosmic rays collide with the air molecules. There is also a 
cosmic neutrino background. Though neutrinos are abundant in the universe, it is an irony 
that they are the least understood. In the Standard Model, neutrinos have no mass, spin 
^ and carry no electric charge. They take part only in weak interactions, which makes 
them extremely difficult to detect as their cross-sections are much smaller compared to 
electro-magnetic cross-sections Q. 

Of the above properties of the neutrinos, the spin and charge of the neutrinos are ex- 
perimentally well known. The mass of the neutrino is relatively unknown. Experiments 
which put kinematic limits on the neutrino mass directly are difficult to conduct and put 
weak limits Q. However, the abundant sources of neutrinos, the stars and the atmosphere 
help us in understanding the properties of neutrinos further. As mentioned earlier, neutrino 
cross-sections are extremely small. Thus to detect them one would need in addition to large 
fluxes (which are naturally provided by these sources), large detectors too. These detec- 
tors measure the number of neutrinos produced for example, in the cosmic ray showers in 
the atmosphere. The results from these experiments over the years consistently pointed to- 
wards the phenomena of 'neutrino oscillations', which in turn indicates existence of neutrino 
masses. 

The above neutrino oscillation experiments can be broadly divided as a) solar neutrino 
experiments, b) atmospheric neutrino experiments and c) laboratory experiments, depend- 
ing on the source of neutrinos. Below, we consider some details of these experiments and 
how they determine the pattern of the neutrino mass matrix. 

2.1 Neutrinos from the Sun and the Solar Neutrino Anomaly 

As we have seen above, neutrino are produced in the stars which burn hydrogen fuel. In 
the Sun, this process produces as a byproduct electron neutrinos, fg- Since the v^, inter- 
act weakly with the solar atmosphere, they can escape the Sun without much changes in 
their flux or energy and thus making it possible to measure their flux and energy on the 
earth. For over thirty years starting from 1967 [^, this flux of Vf. has been measured on the 
earth. It is found that the measured flux is only about one-third of the standard expecta- 
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tion. This discrepancy constitutes the solar neutrino problem. Several other experiments 
Uke Kamiokande, SAGE, GALLEX and super-Kamiokande have since then confirmed the 
existence of this problem, the last one with very high statistics. 

Source Flux Average Neutrino Maximum Neutrino 

(cm~^s~^) Energy (MeV) Energy (MeV) 



pp 


(5.94 ± 0.06) X 10^° 


0.2668 


0.423 ± 0.03 


pep 


(1.39 ±0.01) X 10^ 


1.445 


1.445 




(4.80 ± 0.43) X 10^ 


0.3855 
0.8631 


0.3855 
0.8631 


«B 


(5.15t°:?^) X 10^ 


6.735 ± 0.036 


~ 15 


hep 


2.10 X 10^ 


9.628 


18.778 




(6.051^:^^) X 10^ 


0.7063 


1.1982 ±0.0003 


15q 


(5.32tJ]o^) X 10^ 


0.9964 


1.7317 ±0.0005 


17 p 


(6.331°;?^) X 10^ 


0.9977 


1.7364 ± 0.0003 



Table 2.1 In the above, we present the total flux of the neutrinos coming from various 
interactions along with their average and maximum energies 

In the Sun, the main reactions which produce the energy can be grouped as the pp cycle 
and the CNO cycle. The final outcome of both these cycles can be given as, 

4p±2 e- ^2 z^e+^He2±Q (2.1) 

Here Q = 26.72 MeV represents the energy released in this process. A major part of the 
energy is in the form of photons. A small part (~ 2%) is taken by the neutrinos. Neutrinos 
are produced in various intermediate reactions within the pp and the CNO cycles. In the 
pp cycle, these intermediate reactions are named ( after their reactants ) as pp, pep, Be, 
^B, hep. In the CNO cycle the corresponding reactions involve ^^A^, ^^O/^F. The cross- 
sections of these reactions are determined by the standard model weak interaction physics. 
But the total flux of the neutrinos emanating from these reactions depends on the chemical 
abundances in the solar interior which in turn is dependent on the solar physics. 

These total fluxes are determined by a Standard Solar Model (SSM) developed by J. 
Bahcall and his collaborators The model assumes that a) the Sun is in Hydrostatic 
equilibrium, b) Energy transport happens in Sun by photons or by convective motions, c) 
energy generation is done by nuclear reactions and d) changes in the chemical abundances 
are solely due to nuclear reactions. With these assumptions, the model then predicts the 
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flux of the neutrinos to be observed at various experiments. It sliould be however noted 
that the flux of the neutrinos also changes with the neutrino energy. Infact, neutrinos 
produced from different reactions would have a different energy spectrum. In the Table 
2.1 we present, the total fluxes of the neutrinos from various reactions as predicted by the 
SSM and their corresponding average and maximum energies. 

Experiment Result Theory Theory 



HOMESTAKE 2.56 ± 0.16 ± 0.16 7.7til ^-^^-om 

Gallex 77.5±6.2t^;^ 129tl 0.60 ± 0.07 

Sage 75.4^^;^ U7tl 0.51 ± 0.07 

Kamiokande 2.80 ±0.19 ±0.33 5.15^^;? 0.54 ± 0.07 

super-Kamiokande 2AO±Omtom 5.15l:J:? 0.465 ± 0.005 

Table 2.2: The results of solar neutrino experiments confronted with the corresponding 
theoretical predictions. The results of Homestake, GALLEX and SAGE experiments are 
expressed in terms of event rates in SNU units (1 SNU = 10^^^ events atom^^s^^ ), whereas 
the results of the Kamiokande and Super-Kamiokande are expressed in terms of the 
neutrino flux in the units of 10^cm~'^s~^. The first experimental error is statistical and the 
second is systematic /j^, 



The first experiment to measure the flux of the neutrinos coming from the Sun was the 
chlorine HoMESTAKE experiment which measured about one-third of the standard expecta- 
tion This experiment has a high threshold (~ MeV) and hence was able to measure fiux 
of neutrinos coming from and ^Be only. It reported a deficit of about one-third of the 
SSM prediction. In 1988, the Kamiokande experiment using a real time water cherenkov 
detector could point out the directionality of the incoming neutrinos and confirmed that 
neutrinos are indeed coming from the Sun. This experiment too had a high threshold (~ 6.5 
MeV) and hence was able to see only the ^B neutrino flux. This experiment reported a 
deficit about one-half of the SSM prediction. In 1992, Gallex and Sage with thresholds 
of order ~ 200 KeV were able to measure neutrinos coming from the pp reactions also, con- 
firming that the source of Sun's energy is indeed through thermonuclear fusion reactions. 
They too measured a deficit of about one-half. Recently, Super-Kamiokande experiment 
announced results with very high statistics a deficit of one-half in the Boron neutrino flux to 
which the experiment is sensitive. In Table 2.2 we have listed the latest results from vari- 
ous experiments along with the corresponding SSM predictions. It should be however noted 
that the SSM is also being constantly updated over the years with latest inputs from various 
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experimental results and computational techniques. The theoretical predictions presented 
in Table 2.2 are based on the latest version of the model called BP98 0. An independent 
confirmation of the Standard Solar Model was recently achieved by the helioseismological 
data which confirms the predictions of the SSM for sound velocities in the solar interior to 
a very high accuracy 

In the above, we have seen that the measured solar neutrino flux is about half the 
expectations based on the Standard Solar Model (SSM). One of the possible reasons for this 
deficit could be the expectations themselves which are based on a model. However, it can be 
shown that independent of the SSM, one still faces problem to explain the data from various 
experimental results. Also model independent constraints like luminosity constraint etc, 
which when used in conjunction with the experimental data would lead to scenarios of the 
Sun which are unphysical or severely constrained by other observations like helioseismology. 
For example, using the luminosity constraint one can show that data from Gallex or Sage 
would fit well if all the neutrinos are coming from pp reactions only. Such a scenario which 
suppresses other sources of the neutrinos is severely constrained by helioseismological data. 
Similarly, a simultaneous fit for Super-Kamiokande data and Homestake data would 
leave no room for the flux of neutrinos coming from 'Be. This result is again severely 
constrained by helioseismological data js], |9|. Thus it looks like independent of Standard 
Solar Model, one still has a problem in explaining the data coming from solar neutrino 
experiments. In the present thesis, we believe that the SSM is correct and there exists a 
solar neutrino problem with respect to the expectations based on the SSM. 

2.2 Neutrinos from the Skies and the Atmospheric Neutrino Anomaly 

The other natural source of neutrinos comes from the cosmic rays. Cosmic rays which 
have led to many fundamental discoveries in particle physics early this century, have now 
contributed to our knowledge of neutrinos significantly. Neutrinos are produced when in- 
tergalactic cosmic rays interact with the air molecules in the atmosphere. The production 
mechanism can be summarised as follows: 

(1) . Primary cosmic rays interact with the air molecules producing kaons (K^) and pions 
(vr±). 

(2) . These pions then decay to form a part of the neutrino (i'^^) flux and muons (/x^). 

(3) . Lastly muons decay to give rest of the neutrino (i'^^) flux and the (W) flux. 

The actual flux which reaches the surface of the earth depends on a number of other 
factors like the properties and composition of the primary cosmic rays, modulations due to 
Solar wind and the geomagnetic field cut-off. These fluxes can be calculated within some 
uncertainties. Whereas the modulations due to Solar wind and the earth's geomagnetic 
field can be incorporated within these calculations, the major uncertainties comes from the 
large errors in primary cosmic ray flux measurements However, one can still make 

strong predictions of the neutrino flux which can be summarised as follows : 

(a) Both fg, and Ue, are produced in these showers. 
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(b) . One can also predict with very less uncertainties the ratio of to Ve in these fluxes. 
This is typically of order: 

'^l^i ~ 2 (2.2) 

Ve + 

However, this would vary from place to place due to the geomagnetic cut-off and with energy 
of the neutrinos. 

(c) . Flux of neutrinos with energy greater than 5 GeV is not affected either by geomagnetic 
cut-off or by solar modulation. 

(d) . One can also predict fairly well the angular dependence of the neutrino flux at a given 
experimental site. Taking the direction of the neutrino to be the direction of the incoming 
cosmic ray|^, the angular dependence of the neutrino flux is characterised by two angles, the 
zenith angle {Q) and the azimuth angle ((/>). The zenith angle dependence is caused by the 
increase in the length of the air column (more probability of decay and thus enhanced flux) 
at angles ~ ^. But since this dependence is symmetric over vr — > vr — 0, this dependence is 
characterised only by |cos0| and would not induce an up-down asymmetry. The geomagnetic 
field can induce both zenith and azimuth angle dependence on the neutrino flux. But, as we 
have seen earlier for neutrino energies ~ 5 GeV, this effect is negligible. Thus these effects 
would not be able to induce an up-down asymmetry in the neutrino flux. 

Though neutrinos from the atmosphere have been detected long time ago at experiments at 
Kolar Gold Fields, India and in South Africa, detailed measurements of these fluxes have 
started only about a decade ago. At the energies of the atmospheric neutrinos, the typical 
detection process is through DIS (Deep Inelastic Scattering): 

y +iV ^ /± + X, (2.3) 

where I is the corresponding charged lepton, N is the nucleon and X is the remnant of the 
scattering. The lepton is then detected either through water cherenkov detectors as in 1MB, 
Kamiokande, super-Kamiokande or through iron calorimeters like Baksan and MACRO. As 
we have seen earlier that even though there are large uncertainties in the primary cosmic 
ray fluxes (~ 20%), the predictions for the ratio, eq.( |2.2| ) are free from such uncertainties. 
Thus the experimental results are usually presented in the form of 'ratio of ratios' given as, 

R=S^, (2.4) 
Rmc 

where 

Rexp(MC) = ^ , I (2-5) 
XVe + UeJ exp{MC) 

One would expect this ratio to be one. In the following Table, we present the results 
from various experiments. 



^ This is a good approximation for neutrinos with energy ~ 100 KeV. 
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Kamiokande 
Kamiokande 



0.60i[j:[]^ ± 0.05 sub-GeV 
multi-GeV 

super-Kamiokande 0.652 ± 0.019 ± 0.051 sub-GeV 
super-Kamiokande 0.668 ± 0.034 ± 0.079 multi-GeV 
1MB 0.54 ±0.05 ±0.11 

Soudan-2 0.61 ± 0.15 ± 0.05 



0.57i|]:j]^ ± 0.07 



Table 2.3 In the above we present results from the various atmospheric neutrino experi- 
ments [|, |lll. 

From the above we see that the measured ratios are much different from unity. This con- 
stitutes the atmospheric neutrino anomaly. 



2.3 Neutrino Oscillations 

In the above we have seen that both solar and atmospheric neutrino experiments measure a 
deficit in the neutrino flux. This deficit can be understood in terms of neutrino oscillations, 
requiring neutrinos to have small masses. The main idea is that similar to — oscil- 
lations, the current eigenstates of neutrinos are not its mass eigenstates. This is expressed 
as, 

( Ve\ ( 



U 



(2.6) 



1^2 
V ^3 J 

where are the mass eigenstates and U is the mixing matrix. The neutrinos are produced 
and detected in their current basis, but they propagate in the physical mass basis. Thus 
a neutrino produced can now change its flavour as it propagates a distance L where it 
is detected. Making simplifying assumptions about relativistic nature of the neutrino, a 



simple survival probability formula can be derived using the schrodinger equation |12|. For 
example, for the electron neutrino survival probability, this is given as 0, 



P — 1 



Ulul^Sin^ 



4E 



(2.7) 



where the LHS represents the survival probability. Am?^ represents the mass squared 
difference between i,j neutrino mass eigenvalues and L represents the distance traveled. 
If there are only two generations involved, this formula reduces to, 

^ For intricacies involving the neutrino oscillation formula, please see 



Neutrino Anomalies and Mass Models 



23 



P.,., = 1 - Sin'29Sin' (^^) , (2.8) 

where now the mixing matrix is represented as a rotation by an angle 6. The solar and 
atmospheric neutrino anomalies can be understood in terms of neutrino oscillations, if the 
mass squared differences of the neutrinos and the mixing angles are of the right order. The 
experimental data now constraints regions in the parameter space of the neutrino oscillation 
formula. The deficit in the data of a particular anomaly can be caused either by oscillations 
in two generations or in three generations. Another possibility is that the neutrinos oscillate 
in to a sterile neutrino. 

CHOOZ and simultaneous solutions 

In the present thesis, we consider that there are only three active neutrinos and the 
solar and atmospheric neutrino anomalies have solutions through oscillations among these 
neutrinos. This would generally need a three flavour analysis of the entire solar and atmo- 
spheric neutrino data. However, assuming a hierarchical pattern for the neutrino masses ^ 
an important constraint on the neutrino mixing matrix comes from the CHOOZ experiment 
which simplifies such an analysis |[T5[| . 

The standard two flavour oscillation solutions for the solar and atmospheric neutrino 
data show an hierarchy in Am^. Assuming neutrino masses also have an hierarchy ^ the 
general three flavour oscillation formulae for the solar and the atmospheric neutrino prob- 
lems decouple and reduce to two flavour oscillation formulae as the mass squared differences 
are not linearly independent 0. These are given as, 

P.... = l-^Ul,Ul,S^n'^^^Y2UUl-Ul,) (2.9) 

P.^u^ = l-^U%U%S^n^[^^^Y (2.10) 

The mixing matrix element ?7e3 is constrained by the CHOOZ experiment to be small. 
The CHOOZ experiment which falls in the category of laboratory neutrino experiments 
[p^ , rules out oscillations for the electron neutrino in the mass squared difference range, 
/S.m? ~ 10~^ with a large mixing ~ 1. For hierarchical pattern of neutrino masses, this in 
turn constraints Uez < 0.22. The above formulae are identical to the 2-flavour oscillation 
formula, eq.( |2.S| ), in the limit of small C/gS- In tliis case, we can concentrate only on two 
flavour solutions to the neutrino anomalies. 

The standard solutions for the the solar neutrino anomaly comprise of three 'regions' in 
the oscillation parameter space of Am^ and Siv?26. One of the regions is called the 'Vacuum 



^For a discussion on the patterns of neutrino mass matrix allowed by the experiments, including degen- 
erate spectra please see [Q. 

^Either m^j ~ ^ "^-1^3 or m,^j ^ <C 771.^3. 
^Amli + Amia + Amig = 0. 
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Oscillations' or 'just so'. The Am^ required is very small ~ 10~^^eV^ and the mixing angle 
~ J. The other two solutions require a much larger Am^, typically of O(10~^eV^). These 
two solutions consider matter effects on neutrino propagation whilst the neutrino traverses 
the distance between the core of the Sun and its surface. They allow for matter enhanced 
resonant conversion (MSW mechanism) |jl^ of the neutrinos from electron species to another 
in specific density regions of the Sun. One of the solutions allows for a large mixing angle, 
^ ~ ? - Large Angle MSW. The other region allows for a small mixing ang le 61 ~ 10-3 and 



is called Small Angle MSW |18]. In the case of atmospheric neutrino anomaly, the results 



are much more constrained. The analysis of super-Kamiokande results [19| allow solutions 



for regions in Am^ ~ 10 ^ and Siin?29 ~ 1. These results are summarised in the table 
below: 



Anomaly 


Solution 




Am^ 
eV^ 


Sin^2d 


Solar 


SAMSW 


(4- 


- 12) 10"^ 


(3 - 11) 10-3 




LAMSW 


(0.8 


- 3) 10"5 


0.7-1. 




vo 


(6- 


11) 10-" 


0.47-1. 


Atmosphere 




(4 


- 6) 10-3 


1 



Table 2.4 In the above we present constraints on 2 flavour oscillation solutions for solar 
and atmospheric neutrino anomalies. 

Whereas the above results have been valid during the period where much of the thesis 
work has been done, recent results from the super-Kamiokande experiment do not favour 
some regions of the parameter space presented above. These results have been systematically 
analysed recently by [20| and we summarise them in the following table |^: 



Anomaly 


Solution 


Am^ 


tar?Q 










Solar 


SAMSW 


(1 - 10) X 10"*^ 


(1 - 20) X 10-^ 




LAMSW 


2 X 10-^ - 10-3 


0.2 - 3. 




LOW-QVO 


6 X lO-'^ - 5 X 10-^° 


0.1 - 7. 


Atmosphere 




(1.6 - 5.4) X 10-3 


0.43 - 2.7 


CHOOZ 






;S 0.043 



Table 2.5 Constraints on oscillation parameters in the light of recent super-Kamiokande 
data are presented above [p0|]. 



^These results are presented in terms of tan S as a convenience for the 'dark side' of the parameter space 
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As we have see from the above Table, the vacuum oscihation solution is no longer 
favoured by the latest solar neutrino data. This is because of an independent observable 
called day-night spectrum which aids in distinguishing different possible oscillation scenarios 
for the solar neutrinos. A 'flat' day-night spectrum as such observed would be difficult to 
reconcile with the vacuum oscillation solutions. However, there are other solutions of the 
type 'Quasi- Vacuum Oscillations (QVO)' which contain small matter effects are allowed by 
the recent data. Most of the work in this thesis was done much before the recent results 
from super-Kamiokande were announced. We will comment on the implications due to the 
new data at relevant places. 



2.3.1 Evidence for neutrino oscillations 

Till now we have been assuming neutrino oscillations as a solution to the solar and at- 
mospheric neutrino anomalies. This assumption can be verified in the experiments both 
with laboratory sources as well as natural sources. One of the first indirect experimental 
observation of neutrino oscillations has been reported by the super-Kamiokande experiment 



recently |22]. As we have seen earlier, the incoming atmospheric neutrino flux is indepen- 
dent of the zenith angle, for energies ~ 5 GeV. This is because the geomagnetic field would 
not affect the primary cosmic ray fiuxes at these energies. Zenith angle dependence can 
also be detected if the involved neutrinos are undergoing oscillations. From eq.( |2.7| ), we 
see that oscillation probability varies with the distance traveled by the neutrino L and the 
energy of the neutrino E. Neutrinos coming from the atmosphere to the detector in the 
super-Kamiokande experiment travel distances within a range of 15 Kms - 13,000 Kms. The 
first number is for the neutrinos which are entering the detector vertically (zenith angle ~ 
) whereas the latter number is for neutrinos which travel vertically upwards (zenith angle 
~ vr). If there are no neutrino oscillations, the number of neutrinos entering the detector 
vertically downwards should be equal to the number of neutrinos traveling upwards. This 
is quantified by the up-down asymmetry given as 

A. = E-ZJl (2.11) 

where U represents the number of upward moving neutrinos and D the number of downward 
neutrinos. Super-Kamiokande experiment has measured this number for the case of muon 
neutrinos and found it significantly different from zero. The actual experimental number is 

= -0.296 ±0.048 ±0.01 (2.12) 

The non-zero up-down asymmetry provides first indirect evidence for the existence of neu- 
trino oscillations. 

At present, we do not have a strong experimental signature for the existence of neu- 
trino oscillations for solar neutrino experiments. The major reason being that different 
experiments measure different regions in energy of the solar neutrino spectrum. While 
Super-Kamiokande has had extremely high statistics, it has been sensitive only to high 
energy neutrinos from and hep neutrinos only. Most of the neutrinos (~ 97%) from 
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the Sun are at very low energies {pp neutrinos ) which have not been measured with high 
statistics. While efforts are on for the measurement of such neutrinos, results are awaited 
from the neutrino experiments like SNO and BOREXINO which would help us to point at a 



single oscillation solution to the solar neutrino problem and clear the ambiguities |23]. The 
SNO experiment has a facility to accurately measure the neutral current interactions of the 
neutrinos also. Thus, if the electron neutrinos from the Sun are converted to some other 
type of neutrinos (active), then the flux of the other type of neutrinos can be detected by 
neutral current interactions. This way, we would get a measurement of the total neutrino 
flux and of course a clear signal of neutrino oscillations. The Borexino experiment would 
be able to measure ^Be flux accurately. This would help to conduct a model independent 
analysis of the solar neutrino flux. Thus, in the next few years we hope to have a single 
solution to the solar neutrino problem with accurate values of Am^ and Siv?29. 

LSND and KARMEN 

Till now we have concentrated only on solar and atmospheric neutrino anomalies. In ad- 
dition to these anomalies, there also has been an laboratory experiment called the Liquid 
Scintillator Neutrino Detector (LSND) which has reported evidence for neutrino oscillations. 
One characteristic feature which differentiates the LSND from the experiments discussed so 
far is that it is an 'appearance' experiment. The LSND has reported evidence for i)^ Vf. 



transitions by observing an excess of Vf. events |24] and thus a direct evidence for the ev- 
idence of neutrino oscillations. In terms of oscillation parameters, these results require 
Am^ ~ leV'^, Sin'^29 ~ 10^^. Whereas the above results have been from experiments with 
^ decay at rest, later experiments with /i decay in flight have been consistent with their 
earlier results. Incorporation of the LSND results would require an additional neutrino 
species as it has a characteristically different mass-squared difference which cannot be 
incorporated within three neutrino species which allow only two mass-squared differences. 

The KARMEN experiment was originally built to verify the LSND results. The KAR- 
MEN experiment found negative results for the some of the allowed parameter space of the 
LSND results. However, it was realised that the KARMEN experiment would not be able 
to verify the entire parameter space allowed by the LSND experiment Q. In this thesis 
where we consider only three active neutrino species, we do not consider the results from 
LSND. 



2.4 Mechanisms of Neutrino Mass Generation 

We have seen that neutrino anomalies can be understood in terms of neutrino oscillations 
if the neutrinos have masses 0. In the hierarchical scenario, typically there are two distinct 
mass scales needed for simultaneous solutions of solar and atmospheric neutrino anomalies. 
The atmospheric mass scale \/Aa is ~ 0.1 eV, whereas the solar scale, \/As can be either 
~ 10~^ eV or ~ 10^^ eV. The mixing matrix has atleast one large mixing in the 2-3 sector. 

^See however, Barenboim and Scheck for scheme involving only three neutrinos. 

* Alternative scenarios are also considered in literature [E6|. We do not consider them here. 
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Thus we see that neutrino mass spectrum is characteristicaUy different from the quark mass 
spectrum. 

As we have seen in Chapter 1, the Standard Model has to be extended in some sector 
to incorporate neutrino masses. In addition to explaining the smallness of the neutrino 
masses, these models should also incorporate large mixing as required by the solutions to 
the neutrino anomalies. Some of the extensions may also have experimental signatures which 
may be stringently constrained. Thus a consistent model of neutrino masses should be able 
satisfy all the above requirements. Many such models have been proposed in literature [ |27| . 
Here instead of considering specific models, we consider two of the most popular mechanisms 
which are often used in literature to build models. 



2.4.1 See-Saw Mechanism 

The most natural way of generating neutrino masses in the Standard Model is to add right 
handed neutrinos in to the model and demand neutrinos attain Dirac masses in the same 
spirit as other fermions of the SM. But, in this case, it becomes difficult to justify the 
smallness of the neutrino mass. A natural way of generating small neutrino masses in this 



manner was proposed by Gellman, Ramond, Slansky and Mohapatra, Senjanovic [29| 
and is called the see-saw mechanism. The key idea in to use to the majorana masses for 
the right handed neutrinos, which are naturally allowed by the gauge symmetry to suppress 
the masses for the neutrinos. 

Representing the three left handed fields by a column vector ul and the three right 
handed fields by uji, the Dirac mass terms are given by, 

- C'^ = 9lMdi^r + H.c (2.13) 

where represents the Dirac mass matrix. The majorana masses for the right handed 
neutrinos are given by, 

- = h'kMRUR + H.c (2.14) 
The total mass matrix is given as, 

1 
2 

where the column vector Vp is given as. 



j^total ^ Lt^^Mvp (2.15) 




(2.16) 



The matrix M is given as. 



m4 " (2.17) 
\Md Mr ) 

Diagonalising the above matrix, one sees that the left handed neutrinos attain majorana 
masses of order, 

M'' = -MIM^^Md (2.18) 
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This is called the seesaw mechanism. Choosing for example the Dirac mass of the neutrinos 
to be typically of the order of charged lepton masses or down quark masses, we see that for a 
heavy right handed neutrino mass scale, (left handed) neutrinos masses are suppressed. This 
way the smallness of the neutrino masses can be explained naturally in this mechanism. The 
see-saw mechanism can be naturally incorporated in Grand Unified Theories like SO(IO) and 
also in left-right symmetric models |^^. Though the actual scenarios are model dependent, 
broadly for a range of Mr ~ 10^ — 10^^ GeV, one attains light neutrinos in these models. 
Large mixing |Q and degenerate spectra |Q can also be realised in these models. Recently 
an extensive analysis has been reported which studies proton decay and neutrino masses in 



SO(IO) GUT theories |32]. 



2.4.2 Radiative Mechanisms 

In the above we have introduced additional fermions with a heavy mass scale to generate 
small neutrino masses. One can instead modify the scalar sector of the model, which is 
anyway not well understood. Neutrino masses are now generated radiatively and thus are 
naturally small. This model is called the Zee model after A. Zee who first proposed it [33|. 

Within the Standard Model neutrinos can attain majorana masses by modifying the 
scalar sector. This can be seen by considering the operator, eatLf C L^, where C is the 
charge conjugation matrix and a, h are the SU{2) indices and i,j are the generation indices. 
This can couple to a field transforming either as a singlet or a triplet under SU{2). Models 
with triplet Higgs are considered unattractive as they contribute to p parameter |3^] . Instead 
here we consider models with a singlet field . The coupling of the lepton fields to this 
field is given as, P^eahLf C L)h'^ . 

Though the above coupling violates lepton number one can always conserve it by defining 

to have a lepton number of -2. However if one introduces an additional scalar (j)2 (in- 
addition to the already existing (j)i)^ a new coupling of the form Map(j)a(t>i3h'^ is possible 
which violates lepton number exactly by two units as required for neutrino mass generation. 
In this model which is named as {0i02^} model, neutrino attain masses at the one-loop 
level and thus are naturally suppressed. One interesting fact about this model is that the 
couplings fab are antisymmetric due to the SU{2) metric. This would lead to an interesting 
texture of the neutrino mass matrix whose diagonal elements are zero. Instead of adding an 
additional doublet one can as well add a doubly charged singlet in to the model, k'^^ . In 
this case, neutrinos attain masses at the two-loop level. This is popularly known as Babu 
model in literature [^]. Including both the additional doublet as well as the doubly charged 
singlet would lead to neutrino masses both at the 1-loop level as well as at the 2-loop level. 
Such a scenario may be required to understand neutrino anomalies in these models with 
discrete symmetries like — — Lj- p6[ | . 

In this thesis, we consider an alternative method to generate neutrino masses. In these 
models neutrinos attain masses employing both the 'see-saw type' mechanism as well as 
radiative mechanisms. We will discuss them in detail in further chapters. 
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Chapter 3 



RG scaling and R violation 



Renormalisation Group (RG) Methods have been introduced by Gellman and Low Q to 
tackle the hmitations of perturbation theory in calculating processes at high momenta. 
Since then they have been applied widely in particle physics and in statistical mechanics 
[^. The main idea is that a change in the renormalised IPI function Tr Q due to a change 
in scale of momentum can be compensated by an appropriate change in the value of the 
coupling constant This would require the coupling constants to be scale dependent. 
In a general regularisation process one introduces an arbitrary mass scale fi. Using this 
scale, the coupling constants are now defined on a 'sliding renormalisation scale'. But, 
the unrenormalised function F is still independent of this arbitrary mass parameter @, |5| . 
Requiring the renormalised Tj. to be invariant under a change of scale (in momenta) and the 
unrenormalised F to be invariant under a change of fi would lead to a differential equation 
called the 'Renormalisation Group Equation' This equation now depicts the fact that 
a change in scale can be compensated with a change in the coupling constant. The end 
product of this analysis are differential equations for the coupling constants with respect to 
the momentum scale which are typically of the form: 

= m (3.1) 

where t represents the scale change in momenta. The functions /? are known as the beta 
functions of the theory. Thus, knowing the coupling constant at one particular scale, the 
value of the coupling constant at another scale can be known by integrating the above 
equation. 

Renormalisation group methods are used in field theory to study the asymptotic limits 
of the theory. For example, asymptotically free field theories are those theories for which the 
/3-function vanishes at a high-scale {E — > oo). Various other features of the renormalisation 
group equations (RGE) like fixed points, singularities, sign of the /3-functions would help to 
understand the nature of the theory at high energies. In particle physics, renormalisation 
group studies have been used to mainly probe the nature of physics beyond the Standard 
Model. Used in conjunction with the Grand Unified models, renormalisation group equa- 
tions can probe coupling constant unification at high energies ||6|. The renormalisation 

^ We denote the renormalised function by Tr{E,x, g, y,) where E is the scale of energy of the process, x 
collectively denotes the various momenta, ratios of momenta etc, g are the coupling constants involved and 
/i is the arbitrary scale introduced during regularisation like dimensional regularisation. 
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group studies of the Standard Model have been recently described in In the present 
thesis, we concentrate on Renormalisation Group Studies in the context of supersymmetric 
Standard Model (MSSM) ^ with and without R-parity violation. 

3.1 Renormalisation Group and MSSM 

As we have seen renormalisation group studies allow us to probe the scale of grand unifi- 
cation. Assuming only Standard Model fields to be present up to Mqut, detailed analysis 
have shown that the gauge couplings do not 'meet' at the same point in the momentum 
scale, where the new theory can take over 0]. The MSSM doubles the particle spectrum 
of the SM. This significantly modifies the evolution of the gauge couplings leading to uni- 
fication of the gauge couplings at a scale, Mqut ~ 10^^ GeV. This can be easily verified 
from the renormalisation group equations of the gauge couplings which are presented in 
Appendix for the MSSM case. The solutions of these equations are given as, 

a,(0) = ^ (3.2) 

where the parameters appearing in the above are defined in the Appendix. Choosing ap- 
proximate values for the gauge couplings at Mz as, 

ai(Mz)«^; Q2(Mz)«^; a^{Mz)^^ (3.3) 

A simple algebra leads us to see that 0, 

^«i(0) = a2(0) = a3(0)«^ (3.4) 

for Mgut = 3x 10^^ GeV. The actual analysis, including two loop RGE and threshold effects 
predicts a2,{tz) = 0.129 which is slightly higher than the observed value. The couplings meet 
at the value Mx = 2 x 10^^ GeV The 'exact' unification of the gauge couplings within 
the MSSM may or may not be an accident. But it provides enough reasons to consider 
supersymmetric standard models seriously as it links supersymmetry and grand unification 
in an inseparable manner ||To[| . 

Along with the gauge coupling unification, some grand unified theories also predict 
bottom quark Yukawa Yj, and tau lepton Yukawa, unification at Mqut- Renormalisation 
Group studies of the MSSM Yukawa couplings show — 1^ unification [11 1 for large ranges 



in tan/?. In the Appendix, we have given the RGE for Yt, Yh and 1^- We have not presented 
here the RGE for the first two generation Yukawa couplings. Since the masses of these 
particles are small, the effect of their Yukawa couplings in the RG studies is generally small 
compared to the third generation Yukawa couplings. We have neglected them in most of 
the analyses presented in this thesis except in the chapter 5 (chapter 6) where the second 
(first) generation Yukawa couplings play an important role. 

^The factor | is required for normalisation of the U{1) group. See |ic| ]. 
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Renormalisation Group studies also play an important role in determining the weak scale 
spectra of the MSSM soft sector. As we have seen in chapter 1, supersymmetry is generally 
broken in a hidden sector. The breaking information is then transferred to the visible 
sector through gravitational interactions (mSUGRA). The result of this mechanism is that 
supersymmetry breaking soft terms are added in to the Lagrangian at a high scale ~ Mqut- 
At the high scale 'strong' universality is assumed with the respective soft parameters being 
equal. The weak scale spectrum is modified due to radiative corrections involving Yukawa 
couplings and gauge parameters. The corresponding renormalisation group equations of 
the parameters reflect this phenomena. Using the RGE, one can determine the entire weak 
scale spectrum in terms of the five basic parameters of the model, tan /3, mo, M, A and B. 
Some salient features of these renormalisation group equations are : 

(a) . The gaugino masses evolve in the same manner as the gauge couplings. Thus the 
relation : 

Ml M2 M3 

= = (3.5) 

91 92 93 

holds true for most of the scales up to small two-loop effects [12|. 

(b) . The superpotential parameters are constrained by non-renormalisable theorems and 
thus the RGE of these parameters are proportional to themselves |]l2| . 

(c) . The soft parameters are unconstrained unlike the superpotential parameters. Thus, 
even if the soft parameters are small or zero at the high scale, they can acquire large values 
at the weak scale due to RG scaling. 

Radiative electroweak symmetry breaking: 

As we have seen in Chapter 1, the MSSM weak scale scalar potential requires at-least 
one of the Higgs mass squared to be negative to generate vacuum expectation values to the 
Higgs. Starting with a positive mass squared for the Higgs at the high scale, large radiative 
corrections from the top quark Yukawa can turn the Higgs mass squared negative at the 
weak scale. This can be seen from the solution to the RGE for the rnj^^ which can be 
approximately written as 0, 

<(t^) - <(0) - ^Ml^sY^n (^^) (3.6) 

where the weak scale is characterised by the mass of the Z boson and MgusY characterises 
the typical SUSY breaking mass scale ~ 1 TeV which appears in the equations ( for example, 
squark masses ). This scale together with the large logarithmic factor ~ 66 appears with 
a negative sign in the above solution. Thus even with a positive mass squared at the 
high energies, the Higgs mass squared can turn negative at low energies. This mechanism is 
called radiative electroweak symmetry breaking Within the mSUGRA inspired MSSM, 
with universal boundary conditions, radiative SU{2) x U{\) breaking helps in reducing the 

= Q characterises the high scale in our notation. 
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number of parameters of the theory. Consider the minimisation conditions which we have 
already seen in chapter 1: 



l^' = ''\.^2 o . - o^l (3.7) 



tan2 13 - 1 2' 

Sin2/3 f m|^, + + 2^^ 
= ^-^ (3.8) 

The parameters in'jj^,ni'jj^ are determined at the weak scale by tan/3, mg, M (or equiv- 
alently M2) A. Using the above equations one can determine and Bf^ at the weak scale. 
Only Sign{fi) remains as a free parameter. On the other hand, one can trade to deter- 
mine tan P 01 n to determine mo etc. In the analyses presented in the thesis we have used 
both the sets of parameters. This mechanism is effective for large ranges in tan /3 and other 
MSSM parameters. Radiative electroweak symmetry breaking can also be incorporated in 
models with gauge mediated supersymmetry breaking. In this case, the smallness of the 
logarithmic factor is compensated by the largeness of the squark masses jl4]. 



3.1.1 Analytical and semi-analytical solutions 

In a typical model of supersymmetry breaking, the soft terms are added at a high scale. 
The soft masses and the couplings at the high scale are the input parameters of the model. 
RG evolution would change these masses and couplings at the weak scale significantly. 
However, one would like to express the weak scale soft masses and couplings also in terms 
of the basic input parameters of the model. This would require to have solutions for the 
RGE for the soft parameters analytically. From the RGE given in the appendix we see that 
except for the third generation sfermion masses, rest of the differential equations can be 
solved analytically. The third generation sfermion masses are in general coupled and have 
to be solved numerically. However, Ibanez and Lopez |jl^ have shown that the solutions 
of these equations can be expressed in terms of 'semi-analytical' formulae, thus making 
the dependence on the basic parameters transparent. For example, in mSUGRA inspired 
MSSM with universality at a high scale, the solutions for the sfermion masses would be 
typically of the form: 

mj-(t) = Dmoml + Dm + D^ma mM A + DmA ^ (3.9) 

where mo, M, A, stand for the standard mSUGRA inspired MSSM parameters. The 
coefficients D's are given in terms integrals which can be solved numerically. In this section, 
we present the solutions of the RGE systematically. 

As we have seen earlier, the RG equations for the gauge coupling constants can be solved 
analytically. Since the gauginos also follow similar type of equations, their evolution can 
also be represented with analytical solutions of the RGE. These solutions are given as 



ai{t) = ai{0)zi{t) 
M,{t) = Mi{0)z^{t), 



(3.10) 
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where Zi{t) is defined as: 

Zi{t) ^ 



(1 + 6, di(0) t) 

with i = 1,2,3 over the generations and hi are defined in the Appendix. Similarly, in the 
limit where the first two generation Yukawa couplings are neglected, the equations for the 
corresponding soft masses, eqs.( p.55 ), depend only the gauge couplings and the gauginos. 



The general solutions of these equations can be represented as: 



m)Xt) = m|.(0) + E (^i ^^^(0)') ' (3-11) 



where j = ^Qi,25 f^i,2' -^1,2; -^1,2) -^1,2) ^1 coefficients appearing in the RHS 

of the respective differential equations. The functions ki appearing in the above solutions 
are given as, 

h{t) = ^ (1 - (3-12) 

As an example, the solution for the RGE of the first two generation left-handed squark 
masses is given as, 

16 1 

^Q.a^t) = "^Qi,2(0) + -jMm^h,{t) + 3M2(0)2/C2(t) + -Mi{{)fki{t) (3.13) 

Similarly, the solutions for the rest of the first two generation soft masses can be read 
off from their RGE. Thus we see that the RGE for the first two generation soft masses 
have analytical solutions as the corresponding Yukawa couplings are small and thus can 
be neglected. But as mentioned above, such an approximation is not valid for the third 
generation Yukawa couplings as they are large and can induce large corrections to the soft 
masses through RG evolution. This makes the set of equations coupled and in general has 
to be solved numerically. However, one can still find semi-analytical formulae for the set 



of equations when tan [i is small [15|. In this limit, Yf, and are small compared to the 
Top quark Yukawa and thus can be neglected. Neglecting the Yf, and 1^ couplings in the 
equation for the third generation Yukawa couplings, eq.( 3.47| ), we can arrive at the following 
solution 1 16 1 for : 

where the function Ei{t) is given as, 

Ei{t) = z^{t)-^^/^^^ Z2{t)-^'^^ zi{t)-^^/^^^^ (3.15) 
with the functions Zi{t) defined above. The function Yden{i) is given as, 

Yden{t) = l + QYt{f)) Fi{t) (3.16) 

where ^ 

Fi{t) = j Ei{t')dt' (3.17) 
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lt(0) is the Top quark Yukawa at the high scale which is determined in-terms of the Top 
quark mass at mt [0]. Defining tmt = 2 In {^^^^ , Yt{0) is given as, 

. , YAtmt) 



Yt{tmt) 



Ei{tmt) - 6 Yt{tmt) Fi{tmt) ' 



16 vr^ sinf3'^ 

(3.18) 



In this approximation, the solutions of and Y^ are given as, 

YMMt) 



Yb{t) 



yl/6 
den 



Yr{t) = Yr{0)Es{t) (3.19) 

With ^^(O) and 1^(0) determined in a similar manner as lt(0) in terms of Yb{tmb) and 
YT-{tmtau) with tmb and tmtau having analogous definitions of tmt. The functions E2{t) 
and i?3(t) are given as, 

E,{t) = zi{tf^Ei{t), 

Esit) = Z2{t)-''^''ziit)-^/^'\ (3.20) 

The integral Fi is generally solved numerically. Thus, we call these solutions as semi- 
analytical solutions. In a similar manner, in deriving the solutions for the rest of the 
parameters we completely neglect the effects of bottom and tau-lepton Yukawa couplings. 
This scheme is sometimes known as 'top Yukawa dominance' approximation. In this case, 
we neglect the bottom and tau-lepton Yukawas appearing in the equations for the soft A- 
parameters, eqs. (|3.5C ), /i and -B^ terms, eqs. (|3.53 ), the third generation squark and slepton 



masses, eqs.( p. 601 ) and the Higgs mass terms, eqs.( p.65 ). Since A;, and Aj- are always 



accompanied by their respective Yukawa couplings in these equations, we will not worry 
about their solutions here. The solution for At is given as, 

At{t) = (^MO) - 1^ Ydenit')Kit')dt'^ (3.21) 

where the function k is given as, 

16 13 

K{t) = — d3(0) M3(0)z3(i)' + 3 02(0) M2{0)z2{tf + — di(0) Mi(0)zi(t)2. (3.22) 
3 15 

The integral on the RHS of the above formula is solved numerically. In the case of third 
generation sfermion masses, the observation that only mQ^, rnfj^ and nrj^^ depend on the 
top quark Yukawa, reduces the task significantly. The RGE for m'jj^, and rn'^^ now 
have the same form as those for the first two generation sfermion masses. Thus their 
solutions too have the same form as given in eq. ( 3.11| ) . The coefficients on the RHS of these 
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equations would now form the constants Cf. m^j^ evolves in the same manner as in 
this approximation and thus has the same solution. 

To find the solutions for the remaining third generation sfermion masses, ie, ni^^ {t),mQ^ (t) 
and 'mjj^{t) we observe that though independently they cannot be solved analytically, there 
are combinations of these masses which can be solved analytically. These combinations do 
not contain the top quark Yukawa 

mi{t) = ml^{t)-2ml.^{t) 

2ml{t) = 2 mjj^{t) - 3ml,^{t) (3.23) 

Using the RGE for these three mass parameters, eqs.( tj.60l) , on the RHS of the above, these 
functions are determined as. 



mlit) = mf,3(0)-2m^3(0)-yM|(0)A;3(i)-3M|(0)A;2(0-Y^M2(0)A;i(t) 

mlit) = m^3(0) - 2771^3(0) - 16M|(0)A:3(t) - 9M|(0)fc2(t) - ^Mf{0)kiit) 

16 

2ml{t) = 2 mjj^iO) -Smfj^iO) -16Mi{0)k3{t) - —Mf{0)ki{t) (3.24) 



Using the above relations the differential equations for m'jj,^{t), mQ^{t), mfj^{t) parameters 
can be rewritten in the approximation of top Yukawa dominance as, 

drun^ (t) „ 



dt 

dt 

with the functions Lj are given by, 



dt ' --^^"' -^^^ 



+ 6Yt{t)mtj^{t) = L2{t) 
+ 6Ytit)ml^it) = Lsit) (3.25) 



16 1 

Li{t) = —asit) M^itf + 3a2{t) M2{tf + —ai{t) Mi{tf 
6 15 

- Yt{t) (-^m2(t) + ^mlit) + Mtf^ (3.26) 

L2{t) = ^as{t) M,{tf + ^d,{t) M,{tf 
6 15 

- 2Yt{t)(^-^ml{t) + ml{t) + At{tf^ (3.27) 

Lsit) = 3d2{t) M2{tf + lMt) Mi{tf 

5 

Yt{t) (mlit) - mlit) + Atitf) (3.28) 



- 3 



* Please see Ibanez-Lopez jisj for a difFerent set of combinations which are valid when all the three third 
generation Yukawa couplings are considered. 
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Solving any one of the equations in eqs. (|3.25|), on e can get the solutions of the other 
parameters from the relations given in eqs.( ^.23| , |3.24 ). For example, the solution for the 
equation for m^^^ is given as, 



H2 



m 



H2 



{Q) + j\^{t')Yden{t')dt' 



(3.29) 



The integral appearing on the RHS of the above equation can be solved numerically. 

(t) has 



Expanding the function L^{t) as given in eqs.( 3.26 ), we see that the solution 



m 



Hi 



the desired form of eq.(3.9). Using this solution and the solutions in eqs.(3.24), we can 
derive the solutions for mQ^(t) and mfj,^{t) as per the combinations given in eqs.( |3.2^ ). 
Analogously, we can start with solutions of either itiq^ or mfj^ and arrive at the other two 
solutions using the combinations listed above. 

The parameters and are generally determined at the weak scale. Using the values 
at the weak scale, one finds the values of these parameters at the Mqut- The solutions for 
these parameters at any scale t using weak scale boundary conditions are given as, 



fj.{t) = ix{tz) 



Ydenjt) 
Ydenitz) 



^1/4 



zi{tz) 



3/lObl 



Z2(tz) 



3/262 



(3.30) 



Bf^ takes the form: 



The functions Zi and Irfen have been defined earlier. In the same notation the equation for 

' {Z2{tz) - Z2{t)) 

(3.31) 



B^{t) = B^{tz) - 3 / At{t')Yt{t')dt' - 3- 

Jtz 

3Mi(0) 



62 



5 hi 



{zi{tz) - Zl{t)) 



As mentioned earlier, in deriving the above solutions, we have followed the approach 
of Ibanez and Lopez [Iq]. An alternative approach has been followed by Barbieri et al 



[18|. Recently, these formulae were also given by Carena et al |1£]. All the above solutions 
hold good only in the top-quark Yukawa dominance approximation. This corresponds to 
regions in the tan (j parameter space around 2 — 20 approximately. For larger values of 
tan [5 Yukawa couplings of bottom and the tau leptons would also become comparable and 
subsequently would have to be taken in to account. Recently, Kazakov and collaborators 



have presented analytical solutions to these RGE even in the limit of large tan fi 1 2C ] . 



3.2 R violation and RG evolution 



As we have seen in chapter 1, R-violating schemes in super symmetric theories have a rich 
phenomenology of their own. In these models, additional couplings which violate lepton and 
baryon number are present in the superpotential. The RGE presented earlier for the case of 
MSSM would be significantly modified to take in to consideration these additional couplings. 
The presence of these couplings leads in general to processes which are constrained severely 
by experiments. These experimental limits can be converted into bounds on lepton number 
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and baryon number violating couplings at the low scale. Using the RGE the bounds at the 
low scales can be converted to bounds at the high scale ||2^. Moreover, since these couplings 
also contain flavour violation, the presence of these couplings at the high scale can generate 



additional flavour violating contributions through RG scaling |22, |2^. Recently, infrared 
fixed point studies have also been conducted in the presence of these couplings [p^. 

In this thesis, we concentrate our studies to neutrino mass structure in the presence of 
these couplings. As we will see later, the typical magnitude of the R-violating parameters 
required in this case is very small, for example, ~ 10~^ for the dimensionless A' couplings. 
The presence of such small R-violating parameter would not modify the RG evolution of the 
MSSM parameters significantly |2^. Thus in this thesis, we consider that the RG evolution 
of standard MSSM parameters would not be effected by the presence of R-parity violating 
couplings. 

However the RG evolution of these couplings can significantly affect the neutrino mass 
spectrum in these models. To study the effect of RG evolution on neutrino mass spectrum, 
we divide the R-parity violating couplings as, 

(a) . Bilinear lepton number violation (dimensionful terms) 

(b) . Trilinear lepton number violation (dimensionless ^ijki^ijk terms). 

We study the structure of neutrino mass spectrum in both of these cases separately. As- 
suming only one set of the parameters to be present at a time i.e either only bilinear or 
only trilinear couplings, we derive the RGE for the R-parity violating parameters present 
both in the superpotential and the soft potential. To derive these equations, we have used 



general formulae given by Falck |26|. In this work, the general formulae were derived by 
using effective potential method. Martin and Vaughn presented general formulae up 
to 2-loop order. They used both diagrammatic and effective potential method to arrive at 
their results. We have used their results in chap. 4 where we have derived two loop RGE 
for R-violating parameters. The equations for the R-violating superpotential and soft po- 
tential parameters are presented in the Appendix for each case separately. The equations 



presented here agree with those presented by Carlos and White |22]. We have also checked 
our equations with the equations presented in other papers [^] . 

3.3 RG scaling of dimensionless L-violation and soft sector 

As we have seen earlier, one of the important characteristics of the RG evolution of the soft 
sector is that, even though these parameters are small at the high scale, they can acquire 
large values at the weak scale. This characteristic plays an important role in supersym- 
metry breaking models like Minimal Messenger Model of Gauge Mediated Supersymmetry 
Breaking, where the soft parameters and A vanish at the high scale. However, these 
parameters acquire non-zero values at the weak scale due to RG scaling thus facilitating 
electroweak symmetry breaking This can be clearly seen from the RG equation for 

the parameter, eq.( 3.53| ). In the presence of non-zero gaugino masses and Yukawa cou- 



plings, a non-zero B^ can be generated at the weak scale. Similarly even if the bilinear 
lepton number violating soft parameters are absent at the high scale, they can be generated 
at the weak scale through RG scaling. This would require non-zero lepton number violating 
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trilinear couplings at the high scale |22, 29, 3C]. 

As mentioned earlier, we would divide the R- violation studies interms of dimensionful Cj 
terms and dimensionless A', A couplings. Consider an R- violating model with only dimen- 
sionless trilinear A' or A couplings. In this case, the soft potential at the high scale would 
not contain any bilinear lepton number violating couplings. However, these couplings can 
be generated at the weak scale due to RG scaling. To see this consider the RG evolution 
of the bilinear lepton number violating soft terms. These are represented as B^^ and rri^ff^- 
The corresponding equations for these terms in the presence of only A' or A couplings are 
given in eqs.( |3.69|j3.73 ). From these equations we can see that even if these couplings van- 
ish at the high scale, non-zero values can be acquired at the weak scale in the presence of 
non-zero A' or A couplings. 

The above generation of bilinear soft lepton number violating terms at the weak scale 
has important consequences for neutrino phenomenology. In general, the presence of di- 
mensionless A', A couplings in the superpotential is believed to give rise to neutrino masses 
only at the 1-loop level. But, RG scaling can generate additional contribution to neutrino 
masses which is much larger. This happens as follows. The bilinear soft lepton number 
violating terms generated at the weak scale lead to generation of vevs for the sneutrinos. 
This in-turn leads to a non-zero mixing between neutrinos and neutralinos leading to a 'tree 
level' mass for the neutrino. Thus, RG scaling can induce a 'tree level' mass to the neutrino 
which could be significantly larger than the 1-loop level mass. The effect of this RG scaling 
in a realistic supersymmetric breaking model and its consequences for the neutrino mass 
spectrum have been worked out in chapters 5 and 6. 



3.3.1 Semi-analytical solutions 

In the case of MSSM without R-violation, we have seen that the weak scale soft spectrum 
can be expressed as semi-analytical formulae in terms of the four basic parameters of the 
mSUGRA model, mo, M2, A and tan /3. The coefficients of these parameters involve simple 
integral which can be solved numerically in the limit of low tan f3 where only the top quark 
Yukawa dominates. A similar analysis can be done for the R-parity violating couplings 
discussed so far. As an example, we present here analytical and semi-analytical solutions 
for the case where only trilinear X'-jf^ couplings are present in the superpotential. Moreover, 
we choose there is only one coupling of the form X'-j^ present in the superpotential. The 
analysis can be extended for the presence of additional couplings but would be cumbersome. 

In this limit of neglecting Yj,,!^, the solutions for the A' parameters take the following 
simple form, 

Kjkit) = A^,(0)z3(t)-'/=^''='^2(t)-'/''^^i(t)-'/'°''^ (3.32) 
As mentioned above the soft bilinear lepton number violating parameters, S^. and n^l^u^ 
are generated due to RG scaling in the presence of trilinear A' couplings. The solutions for 
the RGE for these parameters have the generic form of that of the soft masses. The solution 
for the RGE for B^. in this case is given by, 

B,^{t) = Pit) ^B,M + jl (3-33) 
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The functions p{t) and qij{t) are given as, 

q,,{t) = /x(t)A^.(t)/i^,.(t)Qs^(t)+^4.(t)) (3.34) 
For the Tn^_^^ , the RGE leads to the foUowing solution: 

KhA^) = ^Ih,{^) jl'-^3{t')dt' (3.35) 

The function rij{t) is given by, 

Tijit) = ml^{t) + mlit) + 2ml^{t) + 24'^.(t)^j^.(t) + 2ml^{t) (3.36) 

The solution for the trilinear lepton number violating coupling, Afjj which appears in 
the above solutions is given as, 



4;.(i) = 4^.(0) -^|jA,(t')lt(t')'^.3dt'-^Mi(0)ai(0)ei(t) 

16 

- 2M2(0)d2(0)6(t) - yM3(0)Q3(0)6(t) (3-37) 
The functions ^(t) are given as, 

mi = ^,{l-^i(.t)) (3.38) 

Using the solutions for the various soft parameters in the above, one can determine the 
weak scale values of these parameters in terms of the supersymmetry breaking soft parame- 
ters at the high scale which are input parameters of the model. Using these expressions the 
entire neutrino spectrum can now be expressed in terms of the few supersymmetry breaking 
parameters. In all the analyses in this thesis, we have solved these integrals numerically 
using MATHEMATICA© software. 



3.3.2 Numerical Methods 

The above semi-analytical formulae hold true only in the limit of small tan /?. In calculations 
involving models which allow large tan /3, the RGE have to be solved numerically. This 
set of first order differential equations is coupled and non-linear and thus have to be solved 
simultaneously. We have used the standard Rung-Kutta algorithm |^ 32 1 which is 
generally recommended for this purpose. 

■''We have employed fourth order RK method. Sometimes we have also employed standard RK routines 
provided by IMSL. 



RG scaling and R violation 



43 



Yukawa couplings 
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Figure 3.1: Flow chart for numerical analysis of RG equations. 



The methodology employed for the numerical analysis is presented in a pictorial manner 
in Fig.(p2|). Masses of the top quark, bottom quark and the tau-lepton determine the input 
parameters for the third generation Yukawa couplings at the weak scale (characterised by 
the Z boson mass ) for a given tan (3. These will be then used to determine the Yukawas 
at the high scale. At the high scale, the masses and the scalar couplings are given as input 
parameters. Running them down along with the Yukawas determine the weak scale values 
of those parameters. Using the weak scale values of 'm'^_^ and w-^^j one can determine the /x 
and parameters at the weak scale. In turn we use these parameters to determine their 
high scale values. The R-parity violating soft terms contain ^ and in their RGE. The 
values of these parameters are known at Mqut- Running down the entire set of RGE to 
the weak scale one finds the weak scale values of the R-parity violating terms. These are 
then used to find out the neutrino masses. 



3.3.3 Appendix 

Standard MSSM Equations: The MSSM RGE are presented in many papers H]. In 



writing down the below set of RGE we follow |34]. 

Here we use the notation t = 2ln{^^) where Mx corresponds to the high scale and Q 
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corresponds to the low scale. In mSUGRA inspired MSSM Mx is taken to be Mqut = 
3 X 10^^ GeV in our calculations of chapters 5 and 6 whereas Mx varies with A in Gauge 
Mediated Models discussed in chapter 4. tz always corresponds to 2ln{^^), where Mz is 
the Z boson mass. 

Definitions 

We also use the following definitions with gi representing the gauge coupling and h f repre- 
senting the Yukawa coupling of the fermion /. 

= A (3-39) 



167r2 

yf = (3-40) 



Gauge couplings 



dasjt) 

dt 
da2{t) 

dt 



SasW (3.41) 
-02(0 (3-42) 



Gauginos 



daijt) 33, 



= Sa3{t)Ms{t) (3.44) 

^ = -Mm^it) (3.45) 

dMiit) 33 ,,,,,, , 

= --aiit)M^{t) (3.46) 



Yukawa couplings 



"^^'^^^ = Yt{t)(^as{t)+3a2{t) + ^aiit)-mit)-Yi,it)] (3.47) 



dt <.v / y 3 ov / ^ 15 

dYbjt) 
dt 

dYr{t) 



dt 

A-parameters 



= Yi,{t)(^^as{t) + 3a2{t) + ^ai{t)-Yt{t)-GYb{t)-Yr{t)^ (3.48) 
= Yr{t)(^3a2{t) + ^Mt)-^Yb{t)-4Yr{t)j (3.49) 



= -^a,it)M^{t)-3a2{t)M2{t)-^a,it)M^it) 

- 6Yt{t) At{t)-Yb{t) Ah{t) (3.50) 
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= -ya3(i) M^{t) - 3a2{t) Ahit) - ^ai(t) Mi(t) 

- Ytit) At{t) - 6Yb{t) Ab{t) - Yr{t) Ar{t) (3.51) 

= -3a2{t) M2{t) - ^ai{t) M^{t) - 2,Yi,{t) - AYr{t) Ar{t) (3.52) 

)U and parameters 



dn{t) 
dt 


Kt) {-lY,{t) - ^Y.it) - \YAt) + l^t) + ^«i(t)) 


(3.53) 


dB^{t) 
dt 


-3a2(t) M2{t) - ^ai M^{t) -A^Y-,- 3^ n - ZAt Yt 


(3.54) 


First two generation squark and slepton masses 




dt 


16 1 
= -a^{t) Miit) + 3a2{t) Ml{t) + -ai(t) M^{t) 


(3.55) 


dml^^it) 
dt 


= ^dsit) Miit) + ^ai(i) M^it) 


(3.56) 


dt 


= ^as{t)Ml{t) + ^ai{t)M^{t) 


(3.57) 


dml^^it) 
dt 


= Sa2{t) Mi{t) + ^ai{t) M!{t) 


(3.58) 


dm\^^{t) 
dt 


= yai M^{t) 


(3.59) 


Third generation sfermion masses 





dmnAt) 16 o o 1 n 

— 1^ = -a3(t)M|(t)+3d2(t)M|(t) + -ai(t)M2(t) 

- Yt{t) (rnl,^{t) + ml^{t)+ml^{t) + At{tf) 

- Y^{t)(ml^{t)+ml^{t) + ml^{t) + A^{tf) (3.60) 
= ^a3(t)M|(t) + l^«i(t)Mf(t) 

- 2Yt{t) {ml^{t) + ml^{t) + m%^{t) + At{tf) (3.61) 
= ^^a,{t) Miit) + ^^a,{t) Miit) 

- 2n(t) {ml^it) + ml^{t) + m]j^{t) + A(i)') (3.62) 
= 3a2it) Miit) + \a,it) Miit) 

- Y,[ml^it)+m%it) + ml^it)+Alit)) (3.63) 
= -a,it)Mlit)-Yrit){ml^it)+ml^it)+ml^it)+Alit)) (3.64) 
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Higgs mass parameters 



'^"'^^^'^ Sa2(.t)Mi{t) + la,{t)M-tit) 



dt -V / / ■ 5 

- 3n(t) [ml^{t) + ml^{t) + m%^{t) + Ab{tf) 

- Yrit) [ml^it) + ml^{t) + m]j^{t) + Al{t)) (3.65) 

= ^a,{t)Ml{t) + la,{t)Ml{t) 

- Wt{t)[rnl^{t) + rnl^{t) + rnl^{t) + At{tf) (3.66) 
Additional RGE for MSSM with lepton number violation 

Below we present RGE for the lepton number violating superpotential and soft potential 
parameters. These are presented separately for the case of dimensionful and dimensionless 

parameters. In writing them, we have neglected terms higher order in the L-violating 
parameter. This is permissible as we are interested only in small L-violation as required by 
neutrino masses. 

Bilinear R-violation 



^ = 6,(i) Qd2(i) + ^ai(0 - ^F#(i) - ^l*(t)) (3.67) 
"^^'^^^^ = B,M {-l^fiii) - l^tit) + l^t) + ^«i(t)) - ei (3«2(t)M2(t) 

+ ^ai(i)Mi(t) + 3At{t)Ytit) + y,f (i)A^(i)) (3-68) 



dt 



Trilinear R-violation (X') 



dt 



Kjkit) (-Y^fit) - Y/^it) - Yf^ (t) - 2n^(t) - 31^^(i)<5,, + ^asit) 
+ 3a2(i) + - ^ h^umuuimjimk (3.69) 



dt ^ \ 2 ' 2 ' ' 2 ^ ' 10 



i^iHiK'') 2 



j^l^mnitmt) (is,(i) + 4,(i)) (3.70) 
it) {—Ybit) - ^Yrit)^ - ^A^^^(t)/i^^(i) (mj,^{t) 
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+ ml^{t) + 2ml^{t) + 2Al^{t)A^^{t) + 2rn\{t)) (3.71) 

- Aiit)Yi{t) - -Mi{t)ai{t) - 2M2{t)a2{t) - —M^{t)a^{t) (3.72) 
Trilinear R-violation (X ) 



dt 



dBeM 
di 



\ [Xijj{t)Y^^{t)SjsSjk + XjumimsSik) (3.73) 

^^^^ = <HAt) (-^^f w - ^yrit) - Imt)) - K(o 

+ 2 mlit) m%{t) + 2Ayt)Af^{t) + 2m%) (3.75) 

dA^ (t) 7 

= -3Mi(0ai(i)-3M2(t)a2(t)-^Aj^-(i)^f(i)-|4.-(*)^fW (3-76) 
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Chapter 4 



Bilinear R violation and Neutrino Anomalies 



4.1 Introduction: 

Among the Standard Model fermions, neutrinos are the only ones which do not carry any 
electric charge. This leads to the possibility that these particles may be of majorana na- 
ture. Majorana masses to the neutrinos would violate lepton number by two units. To 
generate majorana masses to the neutrinos within the Standard Model, one has to modify 
significantly its scalar sector. On the other hand, Supersymmetric Standard Model natu- 
rally allows for lepton number violation. However, these couplings are generally removed 
by imposing R-parity. To allow the lepton number violating couplings back in to the super- 
potential, one can discard R-parity and impose other symmetries like baryon parity. The 
R-breaking couplings now present in the superpotential would give rise to neutrino masses 
iQ. Another approach is to assume R-parity is spontaneously broken ||2|. This approach 
however leads to a low scale massless mode called majoron Q which can be cured Q. In 
this thesis, we follow the earlier approach where explicit lepton number violating couplings 
are present in the superpotential. 

In the present chapter, we look at a specific model of R-parity violation and study 
the neutrino mass spectrum in these models. As has been noted earlier, lepton number 
violation is characterised either by dimensionful terms or the dimensionless A, A' terms in 
the superpotential. In the present chapter, we consider that only dimensionful ej terms are 
the source of lepton number violation in the superpotential. Such a scenario is sometimes 
known as bilinear R-parity violation in the literature. 

= eiL[H2 (4.1) 

There are several theoretical motivations to consider such a scenario. Theories with 
spontaneous breaking of lepton number can be closely identified with these models if one 
assumes e = hi'R where h is the Yukawa coupling of the gauge singlet and the singlet 
field is represented hy uji Q. The e couplings can be identified as the vevs of the singlet 
fields setting the scale of the lepton number violation. Alternatively one could imagine a 
generalised Peccei-Quinn symmetry whose spontaneous breaking leads to /x and terms 
at the weak scale through dim 5 operators ||5|. Moreover, it is possible to choose the PQ 
charges of the different fields in such a way that the generation of the effective trilinear 
operators is suppressed but bilinear terms are allowed 0. A study of neutrino masses in 
bilinear R violating models is thus theoretically well motivated. Several works have been 
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done to this extent in the recent past ||T|, ^, ^, ^. 

As we have seen in chapter 2, with the advent of super-Kamiokande results on the at- 
mospheric neutrinos with very high statistics, a definitive and clear indication for neutrino 
masses and mixing has been achieved. Though there are still ambiguities associated with 
the solar neutrino experiments, the structure of neutrino mass matrix has achieved more 
clarity with one large mixing as a necessary ingredient if one would like to have simultaneous 
solutions for solar and atmospheric neutrino problems with three active neutrinos. It has 
been known for sometime that neutrino mass matrix and mixing are completely calculable 
in terms of standard supersymmetry breaking parameters and R violating parameters in 
supersymmetry models with R violation. The question then arises is, whether it is possible 
to have a neutrino mass matrix which would be able to provide simultaneous solutions to 
both solar and atmospheric neutrino problems within bilinear R violating models which 
have a minimality of only three R violating couplings. Solutions to neutrino anomalies have 
been extensively studied by Hempfling |^ and Joshipura and Babu ||^] within these models 
though they had not concentrated on simultaneous solutions. These studies have been done 
within the framework of standard supergravity inspired MSSM with universal boundary 
conditions at the high scale ~ Mqut- In the present work we extend the analysis and look 
for simultaneous solutions to solar and atmospheric neutrino problems. We also confine our 
study to models where supersymmetry breaking is communicated through gauge interac- 
tions, in particular to the popular Minimal Messenger Model (MMM) We conduct our 
study 1^] within the framework of Joshipura and Babu [^] which is more analytical. 



4.2 Structure of Neutrino Masses: 



The bilinear lepton number violating couplings have a unique feature associated with them. 
By a simple redefinition of the fields one can remove these dimensionful couplings from 
the superpotential. This is possible because the superfields Li,Hi carry the same quantum 
numbers under SU{2)l x U{\)y and such a redefinition would leave the action invariant 



[11|. For example, considering only one dimensionful coupling in the superpotential €3, a 



suitable redefinition of the type: 



would lead to absence of dimensionful lepton number violation in the superpotential. It 
should be however noted that such a redefinition is not possible in the Standard Model as 
the leptonic fields and the Higgs fields transform as two different representations under the 
Lorentz group. This distinction is removed in the MSSM. 

The above feature of bilinear R violation had led to the general belief that these couplings 
would not be of interest as, by a general redefinition one can transform the dimensionful 
lepton number violation to the dimensionless lepton number violation in the superpotential. 
A crucial feature which works against the above arguments is the presence of supersymmetry 
breaking soft terms. These terms do not preserve the above redefinition in general This 
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complication can completely change the phenomenology of the supersymmetry spectrum 
and neutrino masses in particular. To understand these issues in a better manner we follow 
the framework given by Joshipura and Babu. In this framework most of the features of the 
bilinear R parity violating models are transparent unlike the other studies which were more 
numerical. 

The framework: 



Here we consider a redefinition of the fields at the weak scale. A similar redefinition of the 
fields can be done at the high scale. In models with universal soft terms at the high scale, 
such a redefinition would rotate away the soft bilinear terms too. The superpotential with 
bilinear R-violation is written as : 

W = XiL'iE^H[ + fi'H[H2 + eiL'iH2 + \f^QiD^H[ + \%QiU'rH2 , (4.3) 

where i,j = 1,2,3 represents the generation index. A specific choice of the basis is also made 
such that these fields denote the mass eigenstates of the charged lepton (down quarks) in 
the absence of the terms. We now redefine the leptonic and the Higgs fields in such a way 
that the superpotential ( |4.3D does not contain bilinear ej-dependent terms. The redefinition 
is an orthogonal transformation which we define as follows in the unprimed basis: 



Hi = c^H[ + 53(52(01^2 + siL'i) + C2L3), 

Li = -siL'2 + ciL'i, 

L2 = C2{ciL'2 + SiL[) - S2L'r^, 

L3 = -S3//; + C3(s2(ciL'2 + siL;) + C24); (4.4) 



where 



(4.5) 



ei = ^SiS2S3, £2 = mS2S3, 

£3 = /US3C2, /U' = flC3, 

and fi = (^'^ + ef + €2 + £3)^^^- The consequences of the above redefinition are two fold : 
a) The lepton number violation now reappears as trilinear lepton number violation in the 
superpotential. b)The originally diagonal charged lepton mass matrix Mi now acquires 
non-diagonal parts given by 



Ml 



' A1C1C3 -A2S1C3 
A1S1C2C3 A2C1C2C3 

y A1S1S2 A2C1S2 





-A3S2C3 
A3C2 



(4.6) 



a 



<Hi> 

We denote the mass basis for the charged leptons in the presence of non-zero ej as 



e, fi,T and are defined as: 

Li 



(4.7) 



^Note that we have neglected here a sub-dominant contribution to A4i arising due to sneutrino vevs. 
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where 



OlMiO^ = diagonal. 



(4.8) 



Note that the parameters Aj which denote charged lepton masses when = still need 
to be hierarchical if Mi in (4.6) is to reproduce the charged lepton masses 0. One can 
determine 0] Ol by assuming Ai <C A2 <C A3 and neglecting Ai: 



ol 



where 



Cl — SlA^2 



N2 

^23 



Wi \ 

cos 6'23 - sin ^23 
C2N^^N2^ J \ sin6l23 cos 6*23 / 





-S2C3iVfliV2-l 
r-1 AT-1 



(4.9) 



(1 



s2o2\(l/2) 
■^2*37 ' 



N1C1C2S2 



csrur 



(4.10) 



023 is small even for si^2,3 ~ 0(1). We shall therefore neglect it. The trilinear terms 
generated in the superpotential due to rotation (4.4) assume the following form ||^] in the 
physical mass eigenstate basis of charged leptons: 



W 



tan 63 
<Hi> 



[{Ol)3aLa] {m}pLpel + mfQid\ 



(4.11) 



It should be noted here that the trilinear lepton number violating interactions generated 
by a rotation of bilinear couplings form only a part of the most general set of trilinear 
couplings, especially these couplings are flavour diagonal The above trilinear terms lead 
to neutrino masses at 1-loop [^]. The other contribution to neutrino mass is generated by 
the soft supersymmetry breaking terms in a manner discussed below. 

4.2.1 Sneutrino vevs, neutrino masses and RG scaling 

The soft supersymmetry breaking part of the scalar potential can be written as follows in 
the primed basis: 



(4.12) 



The above soft potential contains bilinear lepton number violating terms which are linear 
in the sneutrino field. The presence of such terms can lead to the generation of vacuum 
expectation values for the sneutrinos. The generation of the sneutrino vev in these models 
should be however contrasted with models where lepton number is broken spontaneously 



^This is because Mi is still diagonal even in the presence of non-zero et 
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[0 . The generation of sneutrino vev in these models is a consequence of the expHcit lepton 
number violating terms present in the superpotential. In the limit these couplings tend to 
zero there are no sneutrino vevs in these models. Hence these theories are not plagued by 
major on problems. 

Once the sneutrinos attain vevs^ neutrinos mix with neutralinos through the kinetic 
terms in the kahler potential leading to a tree level neutrino mass. The neutrino mass so 
generated is proportional to the soft lepton number violating coupling, iJ^. which can be 
~ 0{Msusy)- Thus very large neutrino masses can be generated within these models which 
can be problematic phenomenologically. It has been found that a natural way of suppress- 
ing the neutrino mass would be to incorporate the theory in models where supersymmetry 
is broken at a high scale and the low energy supersymmetry spectrum is determined by 
renormalisation group evolution. As we have seen earlier, in these cases, universality be- 
tween various soft parameters at the high scale is achievable like for example in supergravity 
inspired models. One can extend this universality to the lepton number violating soft terms 
which can lead to a suppressed sneutrino vev at the weak scale. 

The above can be more transparently seen by observing that it is always possible to 



choose a minimum with zero sneutrino vev if B^ = Bi and mfji2 



2 



These equalities 



can be imposed at the high scale through universal conditions. But these are not generally 
satisfied at the weak scale. This is due to the presence of Yukawa couplings which distinguish 
between Higgs and leptons. Thus the RG equations which determine the low scale values 
of these parameters depend on the Yukawa couplings. If one neglects the Yukawa couplings 
of the first two generations then the soft mass parameters related to the first two leptonic 
generations evolve in the same way. Thus, one would have the following non-zero differences 
among low energy parameters: 

,2 \ 



m 



AB 



L,H = iB3-B2,^). (4.13) 
. The latter are obtained by 



These differences would now determine the sneutrino ve 
minimizing the scalar potential expressed in the redefined basis of eqs.(4.4). In this basis, 
one finds: 



'flHl (BHI + C3S3'>3(ASi^ - sIABl) - C2S2S31)2ABl) + H.C. 
-C2S2AmL{s3Hf + C3l'3)D2 + S3C3U3Hf*{AmH - s^AmL) + H.c. 



+ li9'+9"m',\' + \i>i\' + \D2\' + \h\' 



where 



\H'2\r, 



(4.14) 



m 



m 



1^3 
2 

B 



"I" ^sAm// — S2S3AmL, 
^\ — s\AmH — s'^c^AniL, 
m?-!^ + slArriL, 
B^ + sIABh- sIsIABl. 



(4.15) 
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The additional terms do not significantly affect the vev for the standard Higgs since the 
sneutrino vevs and AmH,L, ^Bh,l are suppressed compared to m^susY- This remains true 
even for S3 ~ 0(1). The effect of these terms is to induce vevs for the sneutrinos. These 
can be determined by the stationary value conditions of the above potential as: 



{m$^ + D) 

< ^3 > ~ —^^^^{vii-AmH + slAruL) - ^iV2{-ABH + sIABl)) ■ 
{m$^ + D) 

(4.16) 

We have neglected terms higher order in Ami^H, ABh.l while writing the above equations. 
Note that one of the sneutrino field (= vi) does not acquire a vev in this basis. This vev 
would arise if Yukawa couplings of the first two generations neglected here are turned on. 
It should be noted however that the above sneutrino vevs have been derived from the RG 
improved 'tree level' scalar potential. Minimisation of the full one-loop effective potential 



can significantly modify these values [13|. For the sneutrino vevs, one loop corrections have 
been presented recently by Chun et. a/ [ 14 1 and Hirsch et al [^] . Chun et. al use effective 
potential method whereas, Hirsch et. al use diagrammatic method. In the limit of small 
R violation, as is in our case, Chun et. al have found that the one-loop corrections are 
significant only in the regions of parameter space where the tree level mass is suppressed 
and the loop mass to the neutrinos dominates. We do not consider these corrections here. 
However, as we will see below, the effect of these corrections would be negligible on our 
results as in the present model, we do not encounter such a scenario in the parameter space. 

The sneutrino vevs are zero at the boundary scale Mx corresponding to the universal 
masses. Their weak scale values are determined by solving the relevant RG equations. 
These RG equations can be determined by the general set of RGE presented in chapter 3. 
They are given as. 



b) 



= 3Yb{t){ml,+ml+m% + A, 
= -YAt){ml,+ml + m%+Al) 

dABn 

dt 
dABL 



dt 



3Yb{t)Ab{t) , 

-Yr{t)Ar , (4.17) 



where we follow the same notation as earlier with Yf = pj^, is the mass of the sfermion 
concerned, Aj are the trilinear soft susy breaking terms and t = 2 ln{Mx /Q)- 

The tree level mass matrix generated due to these vev ||T7[| can be written in the 

^In the Appendix we present the details of this derivation 
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physical basis Va as: 



Mo = moOi 



where 



tan 



/o 





\ 





2 






2 


J 


<k = 


< Z^2 > 






< lAj > ■ 





oL 



(4.18) 



(4.19) 



Ol is defined by eqs.(4^) and 



M(cg'+r)(< t>2 + < i>3 >') . 

° 2(-c/iM + 2M2,C/3S^(c + tan2^VK))' ' 
with M2 now representing the weak scale value of the gaugino mass and c = 0.49 

4.2.2 1-loop mass 

The trilinear interactions in eq. ( |4.11 ) lead to diagrams involving squarks and sleptons in the 
loop and generate the neutrino masses at the 1-loop level 1 12]. These contributions depend 
upon the masses as well as mixing between the left and the right handed squarks as well as 
the sleptons. These are however fixed in terms of the basic parameters of supersymmetry 
breaking. In the present case, the trilinear couplings are not independent and are controlled 
by the fermion masses. As a consequence, the dominant contribution arises when the b- 
squark or r slepton are exchanged in the loop. We shall retain only this contribution. 
Let us define: 



61 cos (pb + 62 sin (ph 

62 cos (pb — bi sin (pb 



(4.21) 



Where, 61^2 are the mass eigenstates with masses Mb^^b2 respectively. The mixing angles 
(pr and masses are defined analogously in case of the tau slepton. The exchange of 

b-squark produces the following mass matrix for the neutrinos: 

{Mib)a(3 = "llb(OL)«3(OL)/33- (4.22) 



Due to the antisymmetry of the leptonic couplings in eq. ( [4.1l| ), the exchange of the r 
slepton leads to the following contribution: 



/ o 



M- 



Ir 



V 



il3 
OL13OL23 









D2 





0/ 



(4.23) 



'Ml — cA/2 is used in deriving the above 
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The mixing induced by these contributions is completely fixed by the matrix Ol while the 
overall scale of both these contributions is set by, 

mib^ir = Nc ^3 sin (t>b,T cos (j)b,T In I j^i^ 1 , (4.24) 

where A^c = 3, 1 for the b and f contribution respectively. 

In the above we have considered one loop corrections to be present only through the 
combinations of A A, A' A' couplings. But, to derive the one loop neutrino masses, one has 
to consider the complete 1-loop corrections to the 7x7 neutrino-neutralino mass matrix 
^ . This approach has been followed by ^ and recently by M. Hirsch et.al [15| where the 



complete 1-loop mass matrix is written in the tree level mass basis and re-diagonalised. 
However, the approximation we have made in this work can be justified following the work 
of Chun and collaborators This analysis is based on effective mixing matrix approach. 
Using this method, one can analytically understand the dominant contributions to the 
neutrino sector. The 1-loop corrections to the neutralino sector are of order ~ 6% and thus 
can be neglected. The corrections to the neutrino-neutralino mixing (Dirac-type) part are 
sub-dominant to the corrections induced in the neutrino mass matrix. 1-loop corrections to 
the neutrino mass matrix are through combinations of the couplings X X, X' X' , hr, hr X 



[14|. However the most dominant contribution is through the combination of couplings A A, 
A' A' . 

The total mass matrix including the 1-loop corrections is given by, 

M^ = Mo + Mib + Mir. (4.25) 

We stress that the above Ai^ is in the physical basis with diagonal charged lepton 
masses. This matrix assumes particularly simple form when rotated by the matrix Ol- 



oIMuOl 



/ \ 

nios'^ + miT-N 2^^02320^ moSipc^ + mirN2'^C2S2C^ (4.26) 
y mos^Cif, + mirN2'^C2slcl rrioc^ + mib + mirN2'^S2c'^ J 

This explicitly shows that one of the neutrinos is massless in our approximation of neglecting 
Yukawa couplings of the first two generations. The full mixing matrix analogous to the KM 
matrix is given by, 

U = OlOI. (4.27) 

where Oy is the matrix diagonalising the RHS of eq.( [4.26| ). As we will show the mixing 
angle appearing in 0^ is small due to hierarchy in neutrino masses while, the Ol can contain 
large mixing. Hence, the neutrino masses are determined by the matrix ( [4.26 ) and mixing 
among neutrinos is essentially fixed by eq.(|4.9|). 



^Please see the appendix. 

^Recently an additional diagram has been reported in the hteraturc [Q. We do not consider it here. 
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The above formalism shows that the neutrino masses are greatly suppressed compared 
to the typical supersymmetry breaking scale if AmH,L, ^Bh,l vanish at some scale X. 
The weak scale values of sneutrino vev and hence neutrino masses follow from evolution of 
these parameters. It is clear from eq.( 4.17D that the b and r Yukawa couplings control the 



evolution of sneutrino vev. Similarly, the 1-loop masses following from eq. ( 4.24| ) are also 



controlled by the same couplings. As a result, all the effects of lepton number violating 
parameters ej can be rotated away from the full Lagrangian when the down quark and the 
charged lepton Yukawa couplings vanish. Neutrino masses also vanish in this limit. 

This formalism has been applied to the case of MSSM with universal boundary con- 
ditions in |8|. Before studying this model in the case of gauge mediated models which 
is the work of this thesis, we here review the main results of the works of Hempfling and 
Joshipura-Babu. 

Review of results from MSSM: 

As mentioned earlier, universality at the high scale for the soft parameters leads naturally 
to small neutrino masses at the weak scale. Constrained MSSM i.e, minimal Supergravity 
with universal boundary conditions at Mqut provides an appropriate framework from this 
point of view. The total parameters of the above model in the CMSSM framework are the 
standard CMSSM parameters m, M2, tan/3, A and the sign(//) along with the R-parity 
violating parameters or the three angles si, S2, S3. 

It is found that the loop mass mir contributes negligibly to the neutrino mass spectrum, 
through out the CMSSM parameter space [0, ^. In the limit mi,- <C mi;,, mo approximate 
expressions for the neutrino masses can be given as 

~ mo + mife 

~ \ (4.28) 

mo + mib 

where the approximation <C 1 is used in obtaining the second line of the above. One of 
the neutrinos remains massless in this approximation. The main features of the neutrino 
mass spectrum in this framework can be summarized as follows : 

• For most of the CMSSM parameter space, the tree level contribution dominates over 
the loop contribution giving rise to the hierarchical pattern in the mass spectrum 

• But, there also exist regions in the CMSSM parameter space, where the two contri- 
butions to the sneutrino vevs, eq.( 4.16| ) cancel each other for one particular sign of 



/i. In these regions, the loop contribution, mn, can become comparable to the tree 
level mass, mo and cancellations among mo and mi^ can take place. Here, the two 
neutrinos form a pseudo-Dirac pair with a common mass mQS^ relevant for solutions 
of the neutrino anomalies Q. 
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• The mixing matrix is given by the product OlO^ . Even in the case where Oy allows 
only small mixing, large mixing can be generated from 0_l, as it depends only on the 
ratios of the R-parity violating parameters Q . 

• Even though the neutrino masses are suppressed in these scenarios, they are typically 
of the O(MeV). Thus, the R-parity violating parameters have to be chosen to be 
much suppressed compared to the typical order of the supersymmetry breaking scale 
if neutrino masses are to be of the right order to solve the neutrino anomalies ||^ ^. 

• Numerical results from Hempfling show that solutions for solar (either with MSW 
conversion or vacuum oscillations) and the atmospheric neutrino anomalies can be 
accommodated naturally within these models. 

• Recently an extensive analysis has been reported by M. Hirsch al within the framework 
mSUGRA inspired MSSM. They have reported to have found no solutions for the bi- 
maximal mixing scenarios within these models with universal boundary conditions. 
However, with non- universal boundary conditions one can achieve the required |15|. 



4.3 Gauge mediated models and neutrino masses 

The suppression of neutrino masses due to universality of the supersymmetry breaking 
parameters at some scale not only happens in the MSSM with universal boundary conditions 
(CMSSM), but also in models where supersymmetry breaking is mediated through ordinary 
gauge interactions, which have been introduced in chapter 1. To understand how it happens, 
we consider here in some detail the minimal version of gauge mediated supersymmetry 



breaking [19|. In this case the messenger sector contains only one pair of superfields 
transforming as 5+5 representation of the SU{5) group. They couple to a field S which 
is a gauge singlet. Both the scalar and the auxiliary components of S attain vevs, thus 
introducing a supersymmetric mass scale X = X < S > as well as a SUSY breaking (mass)^ 
differences of order Fs- Models with minimal messenger sector are thus characterized by 
two parameters A = and x = ^. 

All the soft parameters related to MSSM fields are fixed at X in terms of A, x and the 
gauge couplings. The gauginos attain their masses at the 1-loop level whereas the sparticles 
attain their masses at the two- loop level. These masses have the following simple form [20|: 

m^X) = 2A'!^C3aUx) + C2aUx) + h^dj{X)^f{x), 

Mj{X) = a,{X)Ag{x), (4.29) 

where m? represents the scalar masses with i running over all the scalars, whereas, Mj 
represents the gaugino masses with j representing the three gauge couplings. The functions 



f{x) and g(x) have been derived in [21|. Here, 



(4.30) 
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C3 = 4/3,0 for triplets and singlets of S'C/(3)c, C2 = 3/4,0 for doublets and singlets of 
SU{2)l and Y = Q - Ts is the hypercharge. Since the Higgs field Hi and the lepton fields 
Li carry the same quantum numbers, m^^^X) = rn^.^X). The minimal messenger model 
(MMM) is further characterised by the assumption of the vanishing bilinear (B) and trilinear 
(A) soft parameters at scale X. 

The phenomenology of the Minimal Messenger Model has been extensively studied in 
various papers |To| , 22, 2C, 24 1. In this thesis, we follow the work of Borzumati [24|. 
The salient features of this model would not change much even in the presence of the non- 
zero Ei. These can be summarised as follows: Only one parameter essentially determines the 
entire soft spectrum as the dependence of the boundary conditions in eq.( 4.29D on x is very 



mild. In the present analysis we choose x = ^ following [p^ j. The MMM is attractive in view 
of the very restricted structure it offers. But as we will show it turns out to be too restrictive 
if one wants to solve the solar and atmospheric neutrino problems simultaneously. We shall 
thus consider an alternative version on phenomenological grounds in which the boundary 
conditions ( |4.29 ) are still imposed but the value of i?^ at A is not taken to be zero. This 
we call as non-minimal model of Gauge Mediated SUSY breaking. 

Framework for gauge mediated models: 

The neutrino mass framework developed in the earlier section holds good for the Gauge 
Mediated Models too except that the boundary conditions are now defined at the low scale, 
X. The formulae for the neutrino mass matrix remain unchanged. The value of the 
parameter at the weak scale gets fixed through its running. This in turn determines both 
fi as well as tan/3 through the minimisation equations presented in chapter 1: 



m 



Sin2/3 



Hi 



m 



H2 



tan^/? 



tan^ (3 — 
2B^fi 



4 



The presence of induces corrections to these equations, but they are very small as discussed 
below eq.(4.15). The eq.( 4.31| ) therefore holds to a very good approximation. 

In spite of the restricted structure, it is possible to self consistently solve the above 
equations |22, 20, 24| and implement breaking of the SU{2) x U{1) symmetry at low energy. 
Vanishing of the soft i?^ parameter at X makes the analysis of this breaking little more 
involved than in the case of the supergravity induced breaking. One needs to include two 
loop corrections to the evolution of the B^ parameter and also needs to use fully one loop 
corrected effective potential. Details of this analysis are presented in 



, 20 



|. We follow 

the treatment given in |24]. The smallness of the B^ at the weak scale results in this scheme, 
in relatively large value of tan /? and its sign fixes the sign of to be positive. The full 1- 



loop corrected potential was employed in the analysis of |24| but it was found that working 
with RG improved tree level potential also gives similar results provided one evolves soft 



parameters of the supersymmetric partners up to a scale Qq 



(m|(X)m^(X))2. We prefer 



to follow this approach and use the RG improved tree level potential of eq. ([4.14| ) in order to 
determine the low energy parameters at the minimum. We have however included two loop 
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corrections to the RG equations for and AB^^h in determining their values at the 
weak scale. Use of RG improved tree level potential allows us to analytically understand 
the structure of neutrino masses and mixing in a transparent way. 

Neutrino Mass structure: 



The three mass parameters 77io,i6.ir) eqs.( 4.20 , 4.24 ) with appropriate definitions control the 
neutrino masses, mo is determined by solving RG equations ( 4.17] ) along with similar ones 
for parameters occurring in them. We have numerically solved them imposing eq.( [4.29 ) as 
boundary conditions at X. We evolved these equations self consistently up to the scale Qq 
defined above. The mo determined in this manner depends upon fi as well as tan /3, both 
of which are fixed in terms of A. 



Figl 



m„/(GeV s^) 





.xlO« 10= l,2xi[ 



0.2 0.4 0.6 0.8 1 



Fig. 1. The variations of (q^^S) ; a-re shown here with respect to A. mo mildly 

depends upon S2 and the displayed curve is for S2 = 0.8. 
Fig. 2. The function is plotted here with respect to S2- 

The loop contributions are fixed in terms of the squark and slepton masses and mixings 
defined in eq. ( 4.21 ). These are determined from the standard 2x2 matrices involving left 
and right squarks and slepton mixing. The elements in these matrices are also completely 
fixed in terms of A. All the three parameters mo, mif,, mij- depend upon an overall scale S3 
of the R breaking. For small S3 they roughly scale as S3. The ratios mo/s^, mife/s^, mi,-/s3 
are thus determined by A alone |^. We have displayed in Fig. 1 variations of (J^y^-a , and 
with A. One clearly sees hierarchy in the loop and sneutrino vev induced contributions. 
This hierarchy gets reflected in the neutrino masses and one obtains hierarchical neutrino 
masses independent of the overall strength of the R violating parameter S3. The mass ratio 



^mo also depends on S2 very mildly through the vev {U2) which is greatly suppressed compared to (fc's). 
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and hence the hierarchy among neutrino masses are seen to be less sensitive to A. The ttlq 
roughly scales hnearly with A. But since the over all scale of mo is set by S3 which is also 
unknown, a change in A is equivalent to changing S3. Thus, we may use one specific value of 
A and neutrino mass spectrum is then completely fixed by three angles si^2,3 or equivalently 
by the three R violating parameters £«. 



The AmH,L, ABh,l entering mo are determined from the RG equations ( 4.17| ) and are fixed 
in terms of A. For example, 

/i ~ 397.0 GeV , tan /3 ~ 46.39 , 

AniH ~ 192661.23 GeV^ , ArriL 2392.35 GeV^ , 

ABh ~ -14.07 GeV , ABl ~ 0.12 GeV, (4.32) 

when A = 100 TeV. The suppression in Ami, ABl is due to color factors and larger squark 
masses compared to the slepton masses in the model. It follows that the ratio tan (p of the 
sneutrino vev, eq. ( |4.19 ) gets considerably suppressed even when the angle S2 is large. We 



show in Fig 2. the value of s^ as function of S2 for A = 100 TeV. Note that this ratio is 
independent of the values of the other R violating parameters when S3 is small. The small 
value of S(f, leads to very simple expression for neutrino masses. The neutrino mass matrix 
in eq.( 4.26| ) is almost diagonal and one finds: 



^ ~ 4sl . 7- . ~cis2 (7.1 X 10-3-5.6 X 10-3) 



"^^-^^^^ - tan^^^. (4.33) 



"23 



(mo + mib) rriy^ 



The masses are fixed in terms of mo^ib^ir which are determined in terms of A and S3. The 
mass ratio is fixed in terms of S2. The range indicated on the RHS in above equation corre- 
sponds to variation in A from (51 TeV- 150 TeV) and ^23 represents the angle diagonalising 
the matrix in eq.( [4.26 ). 



4.4 Neutrino masses: Phenomenology 

As discussed in the last section, the model considered here implies hierarchical masses and 
large mixing without any fine tuning of the parameters. We now try to see if the predicted 
spectrum can be used to simultaneously reconcile both the solar and the atmospheric neu- 
trino anomalies. The model is quite constrained. Three neutrino masses and three mixing 
angles get completely determined in the model in terms of four parameters namely, A and 
three R violating angles si^2,3- In particular, the angle si characterizing the electron num- 
ber violation does not enter the muon and tau neutrino masses, see eq.( [4.26 ). The mixing 



between neutrinos is largely fixed by the matrix Ol with a small correction coming from 
the angle ^23 ^q. (|0^ ). Thus one has approximately. 



Nl{ciVi - SiC2l^2 + SlS2C3J^3) 
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V^j, ~ A^i(-SiA''2l^l + C1C2N2 + CiS2C3A''2 ^T^i) 

Vr « A^r^iV2"H-S2C3l/2Sl + C2I/3) (4.34) 

Note that s\ (52) determines Vf, — [v^^ — v^) mixing. We must thus require S2 to be large 
in order to account for the atmospheric muon neutrino deficit. The s\ should be small for 
the small angle MSW solution and large for the vacuum oscillation solution to the solar 
neutrino problem. As we now demonstrate these constraints are too tight and one does not 
obtain parameter space in case of the MMM, allowing simultaneous solution for both these 
problems. 




0.2 0.4 0,6 08 1 0.2 0.4 0.6 0.8 1 

s2 



Fig. 3a. The effective v^-v^ mixing angle is plotted here with respect to s\ for A = 
100 TeV in the case of minimal messenger model. Fig. 3b. Contours of Am^ are plotted in 
MMM case, for A = 70 TeV (continuous hues) and A = 150 TeV (dash-dot). For A^, the 
upper (lower) lines correspond to 3 x 10~^^eV^(0.3 x 10^^ eV^). For A5, the upper(lower) 
lines correspond to 12 x 10"^ eV^ (3 x 10"^ eV^). 



4.4.1 MSW and atmospheric neutrino problem in MMM 

The angle s\ can be appropriately chosen to fix the required mixing for the small angle 
MSW conversion. The angle 53 which determines the overall scale of neutrino masses is 
also required to be small. In such a case, the survival probability for the atmospheric 
assumes two generation form and one can take the restrictions on relevant parameters from 
the standard analysis as reported in chapter 2. We have determined the effective v^^ — 
mixing and neutrino masses following from eg. ( [4.27 ) by the procedure outlined in the last 
section. We show this mixing in Fig. (3a). In Fig. (3b), we show the masses for two values 
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of A = 70 TeV, 150 TeV. As seen from Fig. (3a), the S2 = 0.3 — 0.75 leads to the required 
sin^ 26^r = 0.8—1. Fig. (3b) displays the contours corresponding to Ag ~ (3. — 12.) 10~^ eV^ 
and Ayi = (0.3 — 3.) 10"'^ eV^ in the S2 — S3 plane. It is seen that hierarchy among neutrino 
masses obtained in the required region is stronger than needed for a simultaneous solution 
of the solar and atmospheric neutrino problems and there is no overlapping region in the 
■S2 — plane for a combined solution. It is of course possible to solve each of this problem 
separately and get the required amount of mixing as well. 
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Fig. 4a. Contours of Am^ in MMM are plotted for A = 100 TeV. For A^i, the upper 
(lower) line corresponds to 3 x 10~^eV^ (0.3 x 10~^eV^). For A5, the upper (lower) line 
corresponds to 3 x 10"^°eV2(0.5 x 10"^° eV^). Fig. 4b. The effective Sm^ 26*5 and 
Sir?2QA are plotted in the case of the minimal messenger model. The inner lines represent 
contours for 0.9 in both the cases whereas, the outer lines correspond to contours for 0.5 
(0.8) for Sin^2Qs {Sin^29A). 



4.4.2 Vacuum oscillations and atmospheric neutrino problem in MMM 

Unlike in the case of the MSW interpretation, the model can nicely account for the hi- 
erarchies required for solving the solar and atmospheric neutrino problems through vac- 
uum oscillations. This is displayed in Fig. (4a) where we show contours corresponding to 
As = (0.5 - 3) 10"^° eV^ and A^ = (0.3 - 3) 10"^ eV^ in the S2 - S3 plane. Unlike in case 
of the MSW conversion, here there is a large overlap region in S2 — S3 plane which leads 
to the required values for A^^^- Despite this one unfortunately cannot explain both the 
problems simultaneously in a phenomenologically consistent way. This is due to the very 
restricted mixing structure displayed in eqs.( 4.27 ). As discussed in chapter 2, the vacuum 
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oscillation probability in the present case is given by, 

Pe = 1 - ^UlUl^ s\T? (^^) - 2UUl - Ul) , (4.35) 

where the last term comes from the averaged oscillations corresponding to the atmospheric 
neutrino scale. Likewise, the muon neutrino survival probability which determines the 
atmospheric neutrino flux is given by. 



(4.36) 



The amplitude of oscillations is controlled by two effective angles: 

sin2 203 = AU^.U^^ , sin2 29a = "^Ul^U^. (4.37) 



The matrix U appearing above is given by eq.( 4.27| ). Restrictions on these angles required 



for a combined solution of the solar and atmospheric anomaly are worked out in for 
different values of Uf.^. Independent of the chosen values for C/gS one requires, 

sin^ 26*5 = 0.5 - 1 , sin^ 26 a = 0.8 - 1 . (4.38) 

It is possible to choose these angles independently and satisfy above equations in a generic 
three generation case. In our case, the mixings are also determined in terms of si,2 through 
eq. ( [4.9| ) . We have plotted the contours corresponding to restrictions in eq. ( B) in Fig. 4b. 
It is seen that there is no region in si—S2 plane for which the solar and vacuum mixing angles 
can be simultaneously large ruling out the possibility of reconciling atmospheric anomaly with 
vacuum solution in the case of the MMM . 

4.5 Non-minimal model and neutrino anomalies 

We had restricted our analysis so far to the MMM which is characterized by eq.( [4.29| ) 
and the vanishing of the B and A parameters at A. Apart from predictivity, there are no 
strong theoretical arguments in favour of this minimal choice. Such a choice would lead to 
positive B^ at the weak scale which in turn assures fi to be positive as noted above. One 
could consider variations of the MMM which in general result in introduction of additional 
low energy parameters. A class of non-minimal models could contain more complicated 
messenger sector which would influence boundary conditions in eq. ( [4.29| ) . Alternatively, one 
may keep the same messenger sector but introduce some direct coupling between messenger 
and matter fields. This could result in non-zero B^ values at X. In fact, gets generated 



[19, 27] in models which try to understand origin of fi term in gauge mediated scenario 
[28|. B^ may be generated in the absence of messenger-matter coupling if MSSM itself is 
extended. 

Fig. 5. Contours of Am? are plotted in Non- MMM case with ^ < 0, for tan /? = 50 
(continuous lines) and tan /? = 40 (dash-dot) with A = 100 TeV. For A^^, the upper 
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FigS 




(lower) lines correspond to 3 x 10 ^eV^ (0.3 x 10 ^eV^). For A5, the upper (lower) 
lines correspond to 12 x IQ-^ eV^ (3 x 10"^ eV^). 

We shall not consider any specific model here, but would adopt a purely phenomenolog- 
ical attitude to point out possible ways which can allow simultaneous solutions of the solar 
and atmospheric neutrino anomalies. It turns out that prediction of mo is quite sensitive to 
the sign of n term which is fixed to be positive in MMM. This follows from eq.( [4.16| ) which 
shows that two contributions to {v) add or cancel depending on the sign of /i. We may thus 
consider a slightly less restrictive form of the MMM in which we regard the value and sign 
of as independent parameters to be determined phenomenologically. Typically, the B^ 
value remains positive due to the strong Yukawa sector coming from the higher loop terms 
in the RGE for B^. One should thus generate B^ sufficiently negative at the A scale, so 
that, it remains negative even at the weak scale. The typical value of B^ at the A scale 
should be around. 



B^{X) « -12 d3 (Qo) M3(go) Ytt^-1 10-^ M3(Qo) t 



(4.39) 



where, the magnitude of RHS denotes the dominant positive contribution to B^ with 
Yt denoting the top Yukawa coupling. A negative B^ of this order at the X scale can now 
generate a negative at the weak scale. We still assume that the mechanism responsible for 
generation of B parameters does not distinguish between leptonic and the Higgs doublet Hi 
and hence B^ and Bi coincide at the scale X. Due to this, sneutrino vevs are still charac- 
terized by the differences in eq.( [4.13| ) and hence are suppressed. The boundary conditions 
on soft masses are still assumed to be given by eq. (|4.29 ). This particular scenario is now 
characterized by parameters A and B^. As follows from the minimum equation, (4.31), one 



may regard the value of tan (3 and sign of ^ as independent parameters instead of Bn. The 
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magnitude of ^ is determined in terms of these parameters by eg. ( [4.31 ). It is now possible to 



simultaneously account for the atmospheric neutrino deficit and have the MSW conversion 
for the solar neutrinos. This is depicted in Fig. (5) in which we show contours (in S2 — S3 
plane) corresponding to A5 = 3-12 10"^ eV^ = 0.3-3 10"^ eV^ for negative /i and two 
representative values of tan/3 = 40, 50. The magnitude of /i gets fixed by eq. ( [4.31 ) to 379.28 



GeV and 382.03 GeV in the respective cases. The corresponding value at the X scale is 
approximately given as Bf,{X) ^ -13.10 GeV (tan/3 = 40); - 7.06 GeV (tan/3 = 50). It 
is seen that now there is a considerable overlap where two mass scales arise simultaneously. 
As mentioned before, these masses are independent of the value of si which can be chosen 
in the required range namely, 

si = 0.0225 - 0.071 

to allow MSW conversion. The angle 6a relevant for the atmospheric anomaly coincides 
roughly with S2 and as follows from Fig. (5), one can simultaneously account for mixing as 
well as masses needed to solve the atmospheric and the solar neutrino problems. 

4.6 Conclusions 

The supersymmetric standard model contains natural source of lepton number violation and 
hence of neutrino masses. The resulting neutrino mass pattern is quite constrained if source 
of lepton number violation is provided by soft bilinear operators and if the SUSY breaking 
is introduced through gauge mediated interactions. This scenario has the virtue that one 
can obtain hierarchical masses and large mixing in the neutrino sector. The hierarchy in 
masses results from hierarchy in the two different sources of neutrino masses while large 
mixing can be linked to ratio of R violating parameters ej. Overall scale of neutrino mass 
is set by S3 and by the Yukawa couplings of the b and r . Neutrino masses are thus 
naturally suppressed and hierarchical. One however needs to assume relatively suppressed 
R violation, i.e. S3 ~ 10~^ in order to obtain the mass scale relevant for the atmospheric 
neutrino anomaly. This requires that €2 ~ £3 ~ 1 GeV when ~ 1 TeV. 

In case of the minimal gauge mediated model, the three neutrino masses and three 
mixing angles are controlled by five parameters A, x and si^2,3- This proves to be quite 
constraining and does not allow one to obtain simultaneous solution of the solar and atmo- 
spheric neutrino anomalies. In the actual analysis, we chose x = |. A smaller value of x 
would not change this conclusion in view of the mild dependence of the boundary condi- 



tions, eq.( 4.29| ) on x. However, a non-minimal version which allows negative parameter is 



capable of accommodating the MSW effect and atmospheric neutrino anomaly. The number 
of parameters needed are still less than in the models based on the minimal supergravity 
scenario. 

4.6.1 RG equations 



In this subsection we present the two-loop parts of the RGE for the parameters B^, ABh, ABl 
used in the calculations in the MMM model. 
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dBl-^°°P 9 207 o 

= ma2{t)M2{t) + -a2{t)ai{t) {Mi{t) + M2{t)) + —6il{t)Mi{t) 

+ im^{t)M^{t)Yt{t) + ^ai{t)Mi{t)Yt{t) + 16d3(t)M3(t)n(t) 

15 

- l^ai(t)Mi(t)n(t) + ^di(t)Mi(t)y,(t) - 16d3(t)At(t)yt(t) 

- ^a^{t)At{t)Yt{t) - 16d3(t)A(t)n(t) + ^di(t)n(t)A(t) 
15 3U 

- ^di(t)A,(t)n(t) + QA^{t)Y^{t) + 6 (^t(t) + A,{t)) Yt{t)Y^Xt) 

+ lSAt{t)Yt{tf + imAtfMt) (4.40) 



^ = 6n2(t)A^(t) + 3y^(t)y,(t)(A,(t) + A,(t)) 



-di(t)^,(t)n(t) + -ai{t)Mi{t)Yr{t) (4.41) 



2— Zoop 
H 



dt 



3y,(t)n(i) (^r(t) + A(t)) - 18n2(t)A(t) 
3y5(t)yt(t) (^t(t) + Ab{t)) - QYb{t)Ak{t) Qd3(i) - ^ai( 
6y;,(t) (^^d3(t)M3(t) - ^di(t)Mi(t)^ (4.42) 



4.7 Appendix 



We present here the derivation of the tree level mass matrix, eg. ( [4. 18 ). In the presence 



of R- violating couplings neutrinos mix with neutralinos. In the bilinear R-parity violating 
scenarios, this mixing takes place with couplings and the sneutrino vevs. In the Weyl 
basis, the Lagrangian describing the neutrino-neutralino mass matrix is given by, 

C^ass = -]^^IMq^q + H.c (4.43) 
where in the two component notation, is a column vector of neutrinos and neutralinos, 

^0 = [l^e,l^,„l^r,-iXl,-iX3,i^H^,i^H2) (4-44) 

The mass matrix has the following general structure which is of see-saw type: 
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Here the sub-matrix m is of dimension 3x4 and has the following structure: 



m 



( 



V 



31 



V2 



V2 



92 
V2 



LO3 



£2 

£3 ) 



(4.46) 



with the coi representing the sneutrino vevs. M4 is the standard 4x4 neutralino mass matrix 
of the MSSM which has the following form: 



Ma 



( 



V 



Ml 




32 

V2 



V2 





M2 

92 



31 



31 



V2 



\ 



V2 



V2 







(4.47) 



The effective 3x3 neutrino mass matrix is obtained by block diagonalising the above matrix. 
It has the form : 



-mM: 



m 



D 



I Af A1A2 



A1A3 \ 
A1A2 Ki, A2A3 
y A1A3 A2A3 A| ) 



The vector A is defined as, 



A 



^oj — Vie. 



D is given by, 



D = 2{-n Ml M2 + 2 Ml^cpSp{Mi + M2 tan^Ow))- 



(4.48) 



(4.49) 



(4.50) 



Here, the matrix rrie// is written in the basis where ej are not rotated away from the 
superpotential. The rotated form of this matrix is given in the text and takes the form Mq, 



eg. ( 4. 18 ), in the charged lepton mass basis. From nieff we can easily see that it has only 
one eigenvalue, even in the presence of the first generation sneutrino vev. This is a generic 



result of all the R-parity violating models |17, 29 1 
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Chapter 5 



Trilinear R violation and Neutrino Masses 



5.1 Introduction 

In the previous chapter we have been concerned with neutrino masses in the presence of 
only bihnear R-parity violating couplings. As we have seen bilinear R violation provides 
a very constrained framework to accommodate simultaneous solutions for both solar and 
atmospheric neutrino problems. In this chapter we pursue an alternative framework where 
we consider trilinear R parity violation in the superpotential. We will study the structure 
of neutrino mass matrix in this framework and the existence of simultaneous solutions for 
the solar and atmospheric neutrino problems within these models. 

One can immediately see that these models definitely allow more freedom compared to 
the bilinear case as the A' and A are large in number. The A' couplings are 27 in number 
whereas the A couplings are 9 in number. Thus due to the large freedom they offer we would 
expect that simultaneous solutions for solar and atmospheric neutrino problems could be 
accommodated in these models for some definite choice of the parameters. But, in addition 
to bringing in larger freedom these models also bring in large amount of arbitrariness with 
them which would make calculations cumbersome. To overcome this arbitrariness and make 
these models predictive, the following choices are generally made in literature: 

1) . Assume A', A couplings to be hierarchical. This choice is generally made by motivating 
that A' couplings if present in the superpotential are likely to have the same pattern as the 
standard Yukawa couplings which give masses to the fermions Q. 

2) . Assume only a sub-set of the A', A couplings is non-zero. Fix the values of these couplings 
numerically, requiring that the neutrino mass matrix determined by this sub-set of couplings 
would give rise to the required pattern. This analysis can then be extended to all possible 
sub-sets Q. 

In this thesis we do not make any of the above choices. We instead make the assumption 
that all the trilinear couplings are present in the superpotential are free parameters, but they 
are typically of the same magnitude. This is a valid assumption as the R-parity violating 
couplings are required to be suppressed compared to the regular Yukawa couplings to give 
correct order of neutrino masses to understand neutrino anomalies. This assumption has 
been earlier followed by Drees et.al in the same context. In their work, neutrinos attain 
masses at the 1-loop level due to the presence of trilinear interactions. But, as we have 
seen in chapter 3, in any realistic model of supersymmetry breaking, trilinear interactions 
would also modify the weak scale soft potential through RG scaling. This would give 
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additional 'tree level' mass to the neutrinos. In this chapter we reanalyse the neutrino 
mass spectrum in these models taking in to consideration the 'tree level' contribution to 
the neutrino masses also. This would significantly modify the neutrino mass spectrum as 
we will see below. 



5.2 Sneutrino vevs and Neutrino Masses 



For definiteness, we would consider only trilinear A' couplings to be present in the superpo- 
tential. Such an assumption is made for simplicity and we would comment on the inclusion 
of A couplings later on. The superpotential in this case is given as, 

W = h^jQiU;jH2 + hiQdiHi + hlUelHi + ixHiH2 + XijkUQjDl (5.1) 

where we have used the standard notation specified in chapter 1. As we have mentioned 
earlier, the presence of non-zero A' couplings would induce two separate contributions to the 
neutrino mass. One contribution arises due to a soft term linear in the sneutrino vev. This 
soft term gets generated through the loops but leads to a sneutrino vev in the 'tree level' 
of the renormalisation group improved low energy effective potential. This contribution 
indirectly generates neutrino masses through their mixing with the gauginos. A majorana 
mass term for the light neutrinos is also generated directly by the loop diagrams involving 
squarks and sleptons [^, ^. We shall refer to these two contributions as RG induced tree 
level (or simply tree level) and loop level masses respectively. 

The framework 

In the present chapter, we choose to work in the mSUGRA inspired MSSM. In this case 
the soft terms are added at the high scale ~ Mgut- In the limit all A' couplings are chosen 
to be similar in magnitude, the analysis of the neutrino mass matrix would become simpler. 
It is now possible to study the neutrino mass matrix analytically. The two contributions to 
the neutrino mass matrix are given as below. 



5.2.1 RG induced Tree Level Mass: 

In the mSUGRA inspired MSSM framework, the structure of the superpotential dictates 
the structure of the soft SUSY breaking terms. Thus, with only trilinear L-violating inter- 
actions, the soft terms do not contain bilinear terms at a high scale. But as we have seen in 
chapter 3, they are nevertheless generated at the weak scale and should be retained in the 
scalar potential at this scale. The relevant part of the soft scalar potential is now given as: 



Vsoft 



|2 , 2 



-fiB^H^H^-B.MiH^ + h.c 



+ 



(5.2) 



where, we have retained only neutral fields and used standard notation with B^. and 
'^.jj-^ representing the bilinear lepton number violating soft terms. The weak scale value 
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of the lepton number violating soft parameters is determined by the RGE given in chapter 
3. We reproduce them here: 



m 



and the standard MSSM RGE for the parameters on the RHS Q In the above, we have 
confined ourselves with the same notation as in chapter 3. Since we allow only trilinear 
interactions in W^, ''^tiHi ~ = at high scale. In the presence of non-zero A^^j, 
and with the above boundary conditions, the parameters m^t^Hi have low energy 

solutions of the form given in chapter 3. Using the solutions for X'^^j and /ij^-, it is convenient 
to parameterise these solutions as 0, 

R — \' h^K- 
"mliHi = Kpphpn'ip ■ (5.4) 

Here p is summed over generations. The parameters k and k' represent the running of the 
parameters present in the RGE's from the high scale to the weak scale. 

The above soft potential would now give rise to sneutrino vevs, which are described by 
the stationary value conditions for the above soft potential. These are given as, 

< ^» > = 2 ,1 2 \o ' (^-5) 
mj^, + ^"^z cos2p 

where vi and V2 stand for the vevs of the Higgs fields and H2 respectively. The sneutrino 
vevs so generated will now mix the neutrinos with the neutralinos thus giving rise to a tree 
level neutrino mass matrix as has been described in the previous chapter. This matrix has 
the structure: 

2{-cfiM2 + 2MlcpS(^{c + tanel)) ' ^ ' 

where c = and g{g') denotes SU{2){U{1)) gauge coupling. 



^The standard RGE for the soft parameters appearing on the RHS of the above equations do get modified 
in the presence of the A' couplings. But as mentioned in chapter 3, in the limit of very small A' couplings, 
as will be required by our model, the presence of these couplings would not modify the running of the soft 
parameters appreciably. 

The assumption that all A' couplings are similar in magnitude is made at the weak scale. Such an 
assumption at the high scale would in general lead to hierarchical A' at the weak scale. 
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5.2.2 Loop Level Mass 

As we have seen earlier, the trihnear coupHngs in the superpotential would also give rise to 
a loop induced neutrino mass with the down squark and anti-squark pairs being exchanged 
in the loops along with their ordinary partners |^]. This mass can be written as, 

3 

= Y^^ipk^'jkp ^1 sin(j)p cos(l)p Inj^ . (5.7) 

In the above, sin(j)p cos(f)p determines the mixing of the squark-antisquark pairs and M^^ 
and represent the eigenvalues of the standard 2x2 mass matrix of the down squark 
system |^]. The indices p and k are summed over. The mixing sin(j)p cosipp is proportional 
to hp and thus one can write the loop mass as, 

■^ij = Kpk "^'jkp ^k ^p ''^loop ■ (5.8) 

3vi sirKpp cos(l)p M| ^ vf 1 
Exphcitly, m,^^p ^ ^ ~ ^ (5.9) 

with MsusY ~ 1 TeV referring to the typical scale of SUSY breaking. Note that mioop 
defined above is independent of the R violating couplings and is solely determined by the 
parameters of the minimal supersymmetric standard model (MSSM). 

5.3 Neutrino Masses and Mixing 

We now make a simplifying approximation which allows us to discuss neutrino masses 
and mixing analytically. It is seen from the RG eqs.(5.3) that the parameters Kj^, k^^, 



defined in eq. ( ^.4| ) are independent of generation structure in the limit in which generation 
dependence of the scalar masses m^,, mg. and soft parameters ^^^^ and is neglected. 
Since we are assuming the universal boundary conditions, this is true in the leading order 
in which the dependence of the parameters multiplying X'n^j^h^ in eq.( |5.3| ) is neglected. 

dependence in these parameters generated through the gauge couplings will also be 
flavour blind though Yukawa couplings will lead to some generation dependent corrections. 
But their impact on the the conclusions based on the analytic approximation below is not 
expected to be significant ^ The neglect of the generation dependence of Kik, k'-^ allows us 
to rewrite eq.( |5.6[ ) as, 

M-ij = moaiaj , (5.10) 

where a, = X'-j^/^ and mo is now completely determined by the standard MSSM parame- 
ters and the dependence of the R- violating parameters gets factored out as in eq. (|5.8|) . mo 
can be determined by solving the relevant RGE. Roughly, mo is given by, 

rough estimate gives ~ 20% variation in the parameter space due to this approximation. 
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where Mqut = 3 • 10^^ GeV. 

Let us rewrite the loop induced mass matrix as, 



= mioop ai aj + niioop h2 Aij + 0(/if , /if ^) , 



(5.12) 



where Aij = A^23^i32 + -^'i32-^i23 " -^'i22-^i33 - -^i33-^i22- Neglecting 0(/if 2) corrections 
to the loop induced mass matrix, the total mass matrix is given by. 



M 



(mo + niioop) ttittj + mioop ^2* 
« M'ij + mioop h2 Aij . 

The matrix M' can be easily diagonalised using a unitary transformation, 

U^M'U = diag (0, 0, mg) , 

where ms ~ (mo + mioop) ( A3I3 + A2I3 + X^^) 



( 



and U 



C2 S2C3 S2S3 
-S2 C2C3 C2S3 
V -S3 C3 ) 



(5.13) 

(5.14) 
(5.15) 

(5.16) 



with S2 



I „ ^ and S3 = ("i+°2) 



The total mass matrix is now given by, 

U'^M^U ^ ms 

where 

A' = U'^A U and e A',, ; 



diag (0, 0, 1) + eA') 



^loop uD i,D A' 



mi 



oop ' ''2 



(5.17) 



(5.18) 



The last equality follows under the assumption that X!-^ are similar in magnitude and 
mioop < mQ. 

2A' 

If U is a rotation in the 12 plane with an angle 61 defined by tan 29i = , then. 



U'^U^M^UU' 



( e<5i 

ms 6(52 

\ cie^'i3 - sie-4'23 sieA'i3 - cie^'23 
m^ diag{e5i,e52-,'^) ■ 



^22 ^11 
cie^ig - sie^'23 ^ 

1 + e^'33 ; 



(5.19) 



where 5i and 82 arc parameters gcnerically of 0(A'^) if all A^^-^ are assumed to be similar in 
magnitude. The off-diagonal elements will generate additional mixing in the model. But, as 
eA' <C 1, we can neglect these off-diagonal elements. The eigenvalues in this approximation 
are given as, 



m^^ e ms Si ; m^^ ~ e ma (52 



m, 



1^3 



ma 



(5.20) 
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As a consequence, 



nib rno I SjA'l 



(5.21) 



i33, 



The last factor in the above is of 0(1) and the remaining part is controlled completely by 
the standard parameters of the MSSM. Eq.( 5.21| ) may be regarded as a generic prediction 



of the model. It is seen from eqs. (U,|ll|) that typically 



. ~ 10-3 (5.22) 



Thus the neutrino mass ratio in eq.( 5.21| ) is suppressed considerably compared to the hi- 



erarchy obtained when sneutrino vev contribution is completely neglected. As we have 
mentioned earlier, the authors of 1^] have not considered the tree level contribution to the 
neutrinos and thus have achieved a hierarchy of neutrino masses proportional to the Yukawa 
couplings: 

^ = ^, (5.23) 

m^^ nib 

where nis{mb) represents the strange(bottom) quark mass. But, as we have seen in the 
above, the RG induced tree level mass alters the hierarchy drastically as the tree level mass 
is almost thousand times larger than the loop mass for most of the parameter space. The 
exact value of the additional suppression factor due to the tree level mass is dependent on 
MSSM parameters and we will calculate it in the next section. 

The mixing among neutrinos is governed hy K = U U' . The angles appearing in K are 
determined by the ratios of the trilinear couplings and hence can be naturally large. Thus, 
as in supersymmetric model with purely bilinear R violation which we have studied in the 
previous chapter, one gets hierarchical masses and large mixing without fine tuning in any 
parameters. 



5.4 Neutrino Masses : Phenomenology: 

We now discuss the phenomenological implications of neutrino masses, eq. ( pOD and mixing. 
Due to hierarchy in masses, one could simultaneously solve the solar and atmospheric 
problems provided, m^,-^ ~ m^^ ~ 10"^ eV and m^^ ~ 10"^ — 10~^ eV. 

In order to determine these masses exactly, we have numerically integrated the RGE 
eq.(^.3|) along with similar equations for the parameters appearing in them. We have im- 
posed the standard universal boundary condition and required radiative breaking of the 
SU{2) X U{1) symmetry. Solution of the RGE determines both mioop (eq. (|5.8|) ) and mo 
(eq.( ^.10| )). We display these in Fig.l as a function of /i for positive (negative) // and 
tan/? = 2.1, Ma =400 GeV ( JVh =200 GeV ). The ratio is quite sensitive to the sign 
of fi. For /i > 0, this ratio is found to be rather small, typically ~ 10^^ — 10"^, while it can 
be much larger for /x < 0. There exists a region with negative ^ in which > 1. In this 
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region, two contributions to the sneutrino vev in eq. (5.5) cancel and ttiq gets suppressed. 



Barring this region, the is seen to be around ~ 10 ^ — 10 ^ for negative /x leading to 

_ m^mu^ ^ 2 (IQ-^ _ iQ-^) (5.24) 

For m,^^ ~ 10^"*^ — 10^^ eV, one thus obtains m^^ ~ m,^^ ~ 2 (10^"^ — 10"^) eV which is in 
the right range required to solve the solar neutrino problem through vacuum oscillations. 
The typical value of mo ~ GeV found in Fig.l implies through eq.( p.l5D , A' ~ 10"*^. Thus, 
one needs to choose all of this order. Once this is done, one automatically obtains solar 
neutrino scale for some range in the MSSM parameters. 

While, hierarchy needed for the vacuum solution follows more naturally, one could also 
obtain scales relevant to the MSW conversion. This happens for very specific region of 
parameters with negative /li in which two contributions to sneutrino vev, eq. (^]^), cancel. 
As already mentioned, can be in large this region. One then recovers the result of Q, 

namely, eq.( ^.21 ) which allows MSW solution for the solar neutrino problem. 




Fig. la. The absolute values of tree level contribution, mo, the loop level contribution, 
mioop and their ratio "^^"'^ are plotted with respect to fi (positive) for M2 = 200 GeV, A=0 
and tan/3 = 2.1 . The mo and mioop are defined in the text. 
Figure lb. Same as in Fig. (la) but for M2 = 400 GeV and /i (negative). 

We showed in Figures la, lb neutrino mass ratio for specific values of M2 and tan/3. 
Qualitatively similar results follow for other values of these parameters. We have displayed 
in Table. 1 values for the MSSM parameters and what they imply for We have shown 
illustrative values of the parameters which lead to the vacuum as well as MSW solution. 
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The latter arise only for limited parameter range corresponding to cancellations in eq.(p.5|). 
The former is a more generic possibility which arise for larger region with both positive and 
negative values of /i. The MSW solution in the present context will have to be restricted 
to the large angle solution if one does not want to impose any discrete symmetries or 
fine tune A's. It should be noted that this solution is preferred by the present data from 
super-Kamiokande experiments as we have seen chapter 2. 



m 


Ms 




mo 


^loop 


ratio 


GeV 


GeV 


GeV 


GeV 


GeV 


m^2 


1312 


200 


1225 


-16.94 


-0.1689 


1.8 10-^ 


3898 


300 


-3400 


-3.238 


0.0394 


2.3 10"^ 


3921 


350 


-3450 


-2.655 


0.0376 


2.6 lO"'' 


157.2 


300 


-777.5 


-0.0470 


0.0487 


-1.9 10^2 


225.7 


400 


-1038 


-0.0368 


0.0369 


-1.8 10-2 



Tablel: The values of niQ, mioop and ratio of the eigenvalues, for various values of the 
standard MSSM parameters m, M2 and /U for tan (3 = 2.1 , A=0. 

The constraints on mixing matrix K, implied by the experimental results are also easy 
to satisfy keeping all the similar in magnitude. As we have seen in chapter 2, the solar 
and atmospheric neutrino problems can be approximately treated in terms of mixing among 
two generations (12) and (23) respectively as long as Kf,^ is small (as implied by CHOOZ) 
and and m^^ are hierarchical. The constraint on relevant mixing following from the 
experiments can then be translated to constraints on elements of K and are given as follows: 

= S3 sin2 26'2 + 4 sin^^es < 1. 
< 0.18 (5.25) 

= 4(ciC2 - SlS2C3)2 (siC2 -|- 820103)'^ < 1. 

It is possible to satisfy all these constraints by choosing for example, 

C3 = S3 = si = ci = ; S2 = 0.254 

The relative smallness of S2 required here does not imply significant fine tuning and can be 
easily obtained, e.g. by choosing, 

-^133 1 
-^233 ^ 

It is to be emphasized that much smaller value of K^s than at the present level will force 
fine tuning and cannot be accommodated in the model naturally. 



0.6 < 4 Kf,, (1 - Kf,,) 

Ke3 

0.8 < 4 K^, KI2 
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We have so far concentrated on the A'^^^ couphngs alone. The analogous discussion can 
be carried out for Ajjfc couplings appearing in the eq.(^). Here also, the tree level contri- 
bution to neutrino masses will dominate over the loop contribution although the structure 
of mixing matrix will differ slightly due to the anti-symmetry of the couplings Xijk in indices 
i and j. 



5.5 Conclusions 

We have discussed in detail the structure of neutrino masses and mixing in MSSM in the 
presence of trilinear R- violating couplings, specifically Noteworthy feature of the 

present analysis is that it is possible to obtain the required neutrino mass pattern under 
fairly general assumption of all the A^^^ being of equal magnitudes. It is quite interesting 
that hierarchy among neutrino masses is controlled by few parameters in MSSM and is 
largely independent of the trilinear R violating couplings under the assumption that all the 
trilinear couplings are equal in magnitude. Thus one could understand the required neutrino 
mass ratio without being specific about the exact values of large number of the trilinear 
couplings. This 'model-independence' is an attractive feature of the scenario discussed here. 

The key difference of the present work compared to many of the other works is proper 
inclusion of the sneutrino vev contribution. While we had to resort to specific case of the 
minimal supergravity model for calculational purpose, the sneutrino vev contribution would 
arise in any other scheme with A^^^ 7^ at a high scale such as Mqut- Such contribution 
thus cannot be neglected a priori. On the contrary, the inclusion of this contribution makes 
the model more interesting and fairly predictive in spite of the presence of large number 
of unknown couplings. The model prefers simultaneous solutions for solar and atmospheric 
neutrino problems with vacuum oscillations for the solar neutrino problem. However, Large 
Angle MSW also can be accommodated in regions of the parameter space where sneutrino 
vev is suppressed. 
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Chapter 6 



Neutrino Mass constraints on R violation and HERA 
anomaly 

6.1 Introduction 

In the last two chapters we were interested in studying the structure of neutrino mass matrix 
in the presence of exphcit lepton number violation and the possibilities of having simulta- 
neous solutions to solar and atmospheric neutrino anomalies. We have motivated such a 
study by claiming that since lepton number violation naturally occurs in a supersymmetric 
extension of Standard Model it would be the right framework to generate neutrino masses, 
as required for solutions of neutrino anomalies. But lepton number violation has not been 
observed in nature ( except for a possible neutrino mass ), whereas the presence of lepton 
number violating couplings can lead to various interesting processes like tt^uv, the 

violation of weak universality, neutrinoless double beta decay etc. Most of these processes 
have not been observed yet in nature and there exist stringent experimental limits on these 
processes [||. These experimental limits can be converted into limits on lepton number vio- 
lating couplings themselves. The limits on various couplings so obtained have been reviewed 
mil- 

The limits discussed above also have an interesting feature associated with them, which 
is the so called 'single coupling scheme'. In obtaining the limit on a particular coupling, 
one usually assumes that the particular coupling is the most dominant one compared to the 
other lepton number violating couplings This is to facilitate a meaningful analysis. But 
the single coupling scheme has its own drawbacks. The single coupling scheme makes the 
analysis basis dependent [Q]. This basis dependence comes due to the CKM matrix in the 
quark sector. In such a case, one coupling in one basis may correspond to several couplings 
in another basis. 

We have also seen that the presence of trilinear lepton number violating couplings in 
the superpotential would also lead to majorana masses for the neutrinos. The neutrino 
masses are themselves constrained by kinematic limits on them. These limits can be used 
to constrain the lepton number violating couplings present in the superpotential. Thus in 
the present chapter, we do a sort of 'reverse analysis' to what has been done so far in this 
thesis. We choose X'^^ couplings for our analysis and study the limits on these couplings 
from neutrino mass constraints . 

^When we say natural here, we mean it is not disallowed by gauge symmetries as in SM. 

^This is because the electron neutrino mass is most stringently constrained compared to other neutrino 
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As has been mentioned in the earher chapter, neutrinos attain mass both at the tree 
level as well as at the 1-loop level in models with trilinear R violation. The present analysis 
in the literature ||5|, |6| of the constraints arriving from neutrino masses has neglected the 
important and dominant tree level contribution to the neutrino masses which would occur 
when any of the lepton number violating term is present in the superpotential. But as 
we have seen in chapter 5, in any realistic model where supersymmetry is broken at a 
high scale, RG evolution strongly modifies the neutrino mass by giving rise to a tree level 
mass to the neutrino, which we call 'RG induced tree level mass'. Incorporation of this 
additional contribution can change the already obtained ^ limits significantly. The aim 
of this chapter is to systematically derive these constraints and discuss its implications. 
The restriction to only A'^^j^ couplings is also interesting because these couplings have been 
relevant in attempts to understand the HERA anomalous results, which we will discuss 
later in the chapter. 

6.2 Basis choice and definition of : 

Since we consider only trilinear R- violation, the lepton number violating part of the super- 
potential is given as follows: 

= A^fc L'.Q'^D'^ (6.1) 

The above is written in the current basis of the quarks. But, as has been mentioned above, 
in order to meaningfully constrain the trilinear coupling, it is usually assumed that only 
a single coupling is non-zero at a time. While the physics implied by these couplings is 
basis independent, the said assumption makes the constraints on basis dependent since 
a non-zero A' in one basis correspond to several non-zero A''s in the other. The trilinear 
couplings in eq. ( |6.lD can be rewritten |Q] in the quark mass basis as follows: 

= Kjki-^AKij + eiUj)dl, (6.2) 

where the above terms are now in the quark mass basis and K denotes the Kobayashi- 
Maskawa matrix. Even in the mass basis one could choose a different definition for the 
trilinear couplings: 

Kjk = KjiKik (6-3) 

and rewrite (|6.2|) as, 

Wft = K,k{-^^dj + eiK\.n)dl (6.4) 

With the first choice, a single non-zero can lead to tree level flavour violations in the 
neutral sector while this is not so if only one X^^f. {j ^ k) is non-zero. Clearly, there 
are other equivalent definitions of trilinear couplings which are between the two extreme 
cases given in eqs.( |6.2[) and (|6.4| ). The basis dependence can be clearly seen here as the first 
coupling X'i2i is constrained severely by the neutrino mass but the coupling \i2i is not. 



masses. 
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The generation of RG induced tree level mass in the presence of trilinear lepton num- 
ber violating couplings has already been discussed in chapter 5. However, for the sake of 
completeness we repeat it here again. The crucial terms to our analysis are the soft su- 
persymmetry breaking terms in the scalar potential. The lepton number violating parts of 
these terms can be written as: 

Vsoft = D,D]v - vHl + m^^o + c.c + • • • , (6.5) 

where we have explicitly written only those terms which are responsible for generating the 
electron neutrino {y\ = u) mass and AJ^- = Xiij- Note that the presence of only trilinear 
terms in the superpotential generate only the first terms in eq.( |6.5| ) at the high scale in 
conventional supergravity based scenario. The remaining terms are however generated at 
the weak scale and should therefore be retained. In this basis, the general RGE presented 
in Chapter 3 become, 

"^^^ = - {y,^ - «2 - ^ai) - -^Kkh^ + , (6.6) 



+2 At:,A? + 2 mf + 2 m^') (6^7) 



— ~2^kk^k ~ 2 *: ~ 2 
16 7 

- — asMs - SdaMa - — diMi . (6.9) 

3 15 

The rest of the parameters carry the same meaning as in chapter 3. We have kept only 



leading order terms in AJ^^ in writing eqs.( |6.6| ). The other parameters appearing in eqs.(3.6) 
satisfy the standard RG equations to this order in X'^j. The second terms in eqs. (j.E) 
generate non-zero and m^/^^ at Qq. These terms in eq. (|0|) generate a non-zero vev 
which can be determined by minimizing the full scalar potential. These are given by, 

BuV2-mljjVi 
(^) ~ — o -1 — ! (6.10) 

^ ' m? + iM| COS2/3 ' ^ ^ 

is the soft SUSY breaking sneutrino mass. This in turn leads to the following neutrino 
mass: 

( \ fi{cg^ + g'^){< D 

{m^Jtree - 2^_a^lM2 + 2M^Cf,Sf, (c + tan^dw)) ' ^ ' 

where c ~ 5g''^/3g^, g^ and g'^ are gauge couplings and tan/3 = {H2)/{Hi) = V2/V1. It 
follows from the eqs. (6^, |6.10|) that the {jnp^)tree involves the combination (A^^/i|^)2. The 



trilinear interactions in eq.(|6.2|) lead also to the following m^^ at the one- loop level: 



(y \2 j^,j2 
(rriujioop ~ - -J^^k sin cos </>fc Inj^ , (6.12) 
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where we have impHcitly assumed that only one A'^^j^. is non-zero at a time. and M| ^ ^ 
respectively denote the mixing among squarks dk,df, and their masses. The mixing (pk is 
proportional to m^. As a result just like the tree level mass, the {my^)ioop also scales as 
(A^^/i^)^. ruv^ therefore provides a bound on A^^ which can be converted to a bound on 
A'^;^ from eg. (|6^) . The resulting bounds become stronger with increase in tan /? due to the 
fact that B,y,m1^_^ involve = /vcos(3 . For the same reason, bounds display strong 
hierarchy, typically. 



(A 



kkJmax 



m 



m 



It also follows from eq. (|6.3D that, 



Ikk/rnax 



Kik 



(6.13) 



(6.14) 



It is seen from last two equations that the Ax33, (A']^3]^) is constrained most (least) by niy^. 

Fig 1 a ( /i > ) Fig 1 b ( M < ) 
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Figure la. FCNC constraints on A'121 for a) m = 200 GeV and b) m = 50 GeV for 
tanjS = 40. Neutrino mass constraints on A'121 c) m = 50 GeV, tan(3 = 15; d) m = 200 
GeV, tan/5 = 40 and e) m = 50GeV, tanj3 = 40. f) from neutrino less /3/3 decay 

g) X'lii from neutrino mass constraints for m = 50 GeV and tan (3 = 40. Figure lb. Same 
as of Figure la) but with /.i < 



6.3 Bounds on trilinear couplings : 

We now present bounds on the couplings ^ijk obtained from the experimental limit on the 
neutrino mass ruu^ Q. In the present framework, the neutrino obtains a majorana mass 
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and hence stringent constraint on ruu^ would follow from the non-observation of neutrinoless 
f3(3 decay. We shall use somewhat conservative value of rrii,^ < 2.0 eV in numerically 
discussing the bounds on different couplings which lead to rriy^. In order to obtain bound 
on the relevant coupling, we have numerically solved eqs.(|6.6|) along with the other standard 
equations for the parameters appearing on the RHS of eqs. ( |6.6| ) as discussed in Chapter 3. 

We work with R violating version of the MSSM with a standard set of soft supersym- 
metry breaking terms specified at a high scale near Mqut = 3 x 10^^ GeV. We confine 
ourselves to the scenario with radiative breaking of the SU{2) x U{1) symmetry. All pa- 
rameters in eq.( |6.10| ) are specified at a low scale Qo • We have chosen Qq to be Mz 0. The 
parameters By,m'^^^ are assumed to be zero at Mqut ■ Their values at Qq depend upon 
the standard parameters of MSSM which we take to be gaugino (gravitino) mass, M2 (m), 
tan/3 and universal trilinear strength A. By and rn'^^_^ determine {my)tree in terms of these 
parameters and ^^^^{tz), where, tz = 2 ln( ^'^^^ ). The 1-loop contribution also gets fixed 
in terms of these parameters since (j)k and ^ appearing in eq.(12) are determined using 
the standard 2x2 mixing matrix for dt, d'f. system. 
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Table 1. Limits on single A^^^ following from the electron neutrino mass for ^ = 0, 
tan /3 = 30 and > 0. The limits become stronger for larger tan /? and weaker for negative 
/i. The existing limits mentioned in the table are for the relevant squark mass ~ 500 GeV 
and gluino mass M3 = 500 GeV in case of A'^^^^. 

The bounds on different couplings are displayed in Figs. 1 and 2. Apart from being 
dependent on SUSY parameters, these bounds are quite sensitive to the chosen sign of 
the n parameter. This is due to the fact that for one (namely negative) sign of fi, two 
terms appearing in the sneutrino vev, eq.( |6.10 ) cancel while they add for positive sign. 
The suppression in sneutrino vev occurring in the first case weakens the bound on A'. Fig. 

^ Note that change in Qo can alter some of the bounds obtained by minimization of the tree level potential 
and more detailed treatment should include 1-loop corrected potential, see G. Gamberini, G. Ridolfi and F. 
Zwirner, Nucl. Phys. B331 (1990) 331 
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1 displays bounds on couplings A'^^^^ and obtained by demanding niy^ = {m^^)tree + 
{'n^ue)ioop < 2.0 eV. Curves for three representative values of tan/3 and universal gravitino 
mass m are shown. For comparison, we also display in the same figure the existing bounds 
on these couplings. The most stringent bound on X'm is derived from neutrinoless double 
beta decay Q and on X'121 from the process -n^uv Q . These are shown as function 

of M2 in the same figure using MSSM expressions for the relevant squark masses. It is seen 
that the bounds derived here are comparable or better (in case of larger values of M2 ) than 
the already existing ones. 

Fig 2 {^i > 0) 

1 0-1 
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Figure 2. Neutrino mass constraints on for a) tan[3 = 5 and b) tan[3 = 25; on A'^33 
for tanf3 = 5 c) considering only loop contributions and d) loop as well as sneutrino vev 
contributions are shown for m = 50 GeV . 

Fig 2 shows the bounds on A'^33 and A'^32 for > 0. The A'^33 is constrained most and is 
required to be as low as 3 x 10~^ even for tan/3 = 5. For comparison we also show the limit 
on A'^33 which would follow if sneutrino vev is completely neglected as in previous works 

|6[. It is clear that inclusion of this vev drastically alters the bound. We also display limit 



on \'i^2 which was poorly constrained in the previous analysis [10|. The bounds on other 
couplings X'lj^ not explicitly displayed in above figures can be read off from eq.( 6.14 ). These 
constraints are listed in a table for two representative values m = 50 GeV, M2 = 150 GeV 
and m = 50GeV,M2 = 175 GeV and tan/? = 30. These respectively correspond in MSSM 
to mj^ ~ 474, 553 GeV. It is seen that the neutrino mass considerably improves on the 
existing constraints. This table is given to illustrate the importance of the bounds following 
from the neutrino mass. It should however be kept in mind that these bounds are sensitive 
to the choice of the MSSM parameters a fact which becomes evident through figures (1-2). 
Thus, we have investigated the restrictions on the R parity violating trilinear couplings, 
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specifically the A'^^, implied by neutrino mass. These restrictions have been discussed earlier 
iP, |6| but earlier discussions either neglected very important and dominant contribution due 
to the induced sneutrino vev or did not incorporate the effect of the quark mixing. We have 
attempted here a complete analysis which incorporates these and also pays careful attention 
to the dependence of the bounds on the choice of basis. Unlike the earlier works, we have 
also numerically discussed the bounds in terms of parameters of the MSSM. 

Many of the discussions and in particular the generation of sneutrino vev in the presence 
of a single trilinear coupling would be true in a more general set up, e.g. in case when soft 
terms arise from the gauge mediated supersymmetry breaking. The main point following 
from our analysis is that the electron neutrino mass provides much stronger constraints on 
R violation and associated phenomenology than has been hitherto realized. 



6.4 HERA anomalies 



The HERA experiments at DESY, collide electrons (or positrons) with protons. The elec- 
trons typically carry an energy of ~ 30 GeV, whereas the protons carry an energy of around 
~ 820(920) GeV. In addition to measuring the structure functions of the protons through 
the DIS (Deep Inelastic Process), these experiments are also ideal to verify R- violating inter- 
actions. The experiment consists of two general purpose detectors, HI and ZEUS. In 1997, 
anomalous events have been reported by the HI and the ZEUS detectors |11| at HERA in 
the deep inelastic e^p scattering. They claimed to have found a resonance for a leptoquark 
of mass ~ 200 GeV. This has generated considerable excitement within the community [|l2[ 
as, if substantiated these results could form strong evidence for physics beyond the SM 
and may be for supersymmetry. The work presented here is based on ref.|jl^ which was 
motivated from the results discussed above. However, this evidence for the leptoquark did 
not stand the test of time and the subsequent data from HERA [ 14 1 were consistent with 
the Standard Model expectations. We describe here the work in ref.[^] for the sake of 
completeness, though it is now only of academic interest. To this extent, we consider here 
data reported by HERA in 1997. 

The available data when taken seriously allowed for two possible interpretations: {i) 
The presence of some lepton number violating contact interaction p5l or {ii) production 
of a resonance in the e'^q channel. MSSM with R-violation |12, 16, ^] provides a natural 
theoretical framework to incorporate the second possibility although an alternative in terms 
of a scalar leptoquark |18l is open. 

The super symmetric interpretation of the HERA events assumes that the excess events 
seen at HERA are due to resonant production and subsequent decay of the squark to e'^q. 
Three possibilities have been considered in this context[^, |l^, 17]: e^d/j — > CL,e~^dR 
ii, e.~^SR ^ ii. In analyzing these scenarios 16, 17] it has been implicitly assumed that 
the squark masses are free parameters of the model. While this would be true in the most 
general situation, specific model dependence can alter some of the conclusions. Our aim 
is to show that the very minimal model dependent assumption on the charm squark mass 
necessarily requires large X'121 to understand HERA events and this large coupling by itself 
is ruled out from other constraints like neutrino mass bounds which we have obtained in 
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the previous section. 

The specific assumption that we make and which leads to the above conclusion is that 
the charm squark mass squared is positive at the unification scale. This assumption is true 
in the radiative electro-weak breaking scenario with universal boundary conditions at the 
GUT scale, but it can also be true in a much more general context. We shall first assume 
that the gaugino masses are unified at Mqut but demonstrate later that the removal of 
this assumption does not significantly change the basic conclusion. The argument leading 
to the above conclusion is largely insensitive to the details of the radiative SU{2) x U{1) 
breaking in the MSSM and runs as follows. 

Consider the following R violating couplings: 

= Kjki-'^idiKij + eiUj)dl (6.15) 

The above terms are defined in the quark mass basis and K denotes the Kobayashi-Maskawa 
matrix. The charm squark interpretation of the HERA anomaly requires to be non- 
zero. The HERA data can be explained provided 

0.025 - 0.034 

^121 ~ ^172 (6-16) 

The number in the numerator of eq.(2) is indicative of the required range and depends upon 
the weightage given to the different experiments as well as on the next to leading order QCD 
corrections In the following, we shall take the value 0.025 for the numerator in the 
RHS of eq.(2). B refers to the branching ratio for the squark decay to qe~^. This decay would 
take place through the coupling in eq.(l) itself. B is also influenced by the R conserving 
decays of the charm squark to an s (c) and a chargino (neutralino) . The the 
parameters ^,M2,tan/3 determine B in the MSSM. HERA data can be reconciled if for a 
region in these parameters (i) eq.(2) is satisfied, (ii) A'^21 is consistent with other constraints 



[y, |g, ^ due to R breaking and (Hi) charm squark has a mass around 180- 220 GeV. 

In supergravity based models, the charm squark mass at the weak scale is governed by 
the gauge couplings and the gaugino masses. Its value at Qo = 200 GeV is given in the 
limit of neglecting Kobayashi-Maskawa and cl — cr mixing by ^ 

m?^(Qo) ~ ml^iMcuT) + 8.83M| + 1/2 M| cos 2^ (1 - 4/3 sin^ Ow) (6.17) 

where we have assumed that the gauge couplings and the gaugino masses are unified at the 
GUT scale, Mqut = 3 x 10^^ GeV and chosen as{Mz) = 0.12. The A'h in eq.(3) is the 
value of the wino mass at the weak scale identified here with Mz- The last term in the 



^ The lower limit is obtained when HI and ZEUS data from 1997 run are also included while the upper 
limit corresponds to inclusion of HI data alone. In both cases, 30% increase in the relevant cross section 
due to next to leading order corrections |l9| is assumed. 

^ We also neglect the effect of additional trilinear R violating couplings on the running of the charm 
squark mass. Their inclusion does not significantly alter the charm squark mass even when they are large, 
see e.g. K. Cheung, D. Dicus and B. Dutta [ pl[ . 
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above equation is a (-ve ) contribution from the D-term. It follows that the charm squark 
mass provides strong upper bound on M2 as long as m?^ {Mqut) > 0: 



M2 < 74.04 GeV 



V22O GeV. 



1-Omcos2fi[—- j j 



The branching ratio B is determined in the MSSM by the following widths [22 
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(6.19) 



(I Fl\'^ + \ Fr |2)(m? - - m^o) - Am^m^oReiFRFl) (6.20) 



A2 (m?^,m^,mj+) 
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We have adopted the same notation as in [22|. From the expression for B in terms of 



the above decay widths, and the HERA constraint, eq.(2), one can solve for the allowed 
\'i2i- The contours in the /i — M2 plane for different values of A^2i ^-^^ displayed in fig la 
(tan/?=l) and fig lb (tan/3 = 40). The horizontal lines in these figures show the upper 
bound on M2, eq.(4). We also display, the curves corresponding to two representative values 
of the chargino masses namely 45 and 85 GeV. The later is the present experimental bound 
obtained assuming R conservation. This need not hold in the presence of R violation. It 
is seen from fig. lb that for chargino mass around 85 GeV, the bound on M2 by itself rules 
out charm squark interpretation for large tan [3 independent of the value of A'^21 0- But 
irrespective of the value of tan/3 and the chargino mass one needs very large \'i2i > 0.13 
in order to satisfy the bound on M2 coming from the charm squark mass. This strong 
bound on X'121 arises because of the following reason. For M2 < 74 GeV, at least one of 
the charginos is sufficiently light and contributes dominantly to the cl decay. This reduces 
B and results in large value for A'^21 due to eq.(2). In contrast, the chargino decay is 

^Specifically, the upper bound on M2 can be reconciled with the chargino mass of 85 GeV or more only 
if tan /3 < 2.5. 

^ Reduction in the branching ratio for charm squark decay in case of the minimal model was also noticed 
in IitII. 
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suppressed kinematically for tan/3 ~ 1 if M2 > 200 GeV. This results in smaller allowed 
value as seen from the figure. But these are in conflict with the charm squark mass. 

Let us now see if one could make large consistent with other constraints. The 
strong constraints come from atomic parity violation |^^, the decay K'^ nuv [Q] and the 
electron neutrino mass which we have derived above. The data from Cs on the relevant 
weak charge have been argued 

A'i2i < 0.074 (6.22) 

at 2fT level in conflict with the large value required here. In principle, the extra contribution 
due to charm squark to atomic parity violation can be canceled by a similar contribution 
from the scalar bottom or strange squark but the existing constraints on the relevant cou- 
plings make this cancelation difficult ||2^. Thus, one cannot easily avoid the atomic parity 
violation constraint strictly in the MSSM but this can be done by postulating new physics, 
e.g. the presence of an extra Z ||2^ . 



12 



to imply 




Fig 1 a: The contours (continuous lines ) of constant A121 obtained by imposing HERA 
constraint, eq.(2). The contours are for values 0.05, 0.08, 0.12, and 0.13. The horizontal 
dashed line represents the bound on M2 coming from requiring mg^ = 220 GeV. The vertical 
dash-dot lines represent the bounds on the chargino mass, the upper one for a mass of 85 
GeVand the lower one for a mass of 45 GeV. All the above are computed for tan/3 = 1. 
Fig 1 b: Same as fig la but for tan/3 = 40 and A'i2i=0.05, 0.08, 0.12, and 0.135. 

The other significant constraint comes from the decay -kvv which implies Q 
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The electron neutrino mass also gives similar constraint in the same parameter range. The 
question of choice of the basis becomes relevant in the discussion of these constraints. This 
is particularly so when one assumes only one X'^^^ to be non-zero. For \'i2i defined in the 
mass basis as in eq.(l) the above constraint is unavoidable if rest of the couplings are zero. 
This basis choice is natural from the point of view of interpreting HERA results but is not 
unique. One may redefine the couplings as 

and rewrite eq.(l) as follows: 

Wr = X'.jki-vrdj + eiKlui)dl (6.24) 

HERA result would now require A'121 to be large. If this is the only non-zero X'ij^ then 
there will not be any constraint on A'121 from the neutrino mass or from the — > t:vi> 
decay ||. But eq.(lO) will now generate a contribution to the neutrinoless double beta 
decay which is also severely constrained. Specifically, one has HI 

This clearly does not allow A'121 of 0(0. 1). Thus, notwithstanding basis dependence one 
has problem in accommodating large value for the relevant coupling. An alternative is to 
allow more than one non-zero X'ijk- It is seen from eq.(lO) that A'lji (j=l,2,3) contribute 
to the neutrinoless (5(5 decay and simultaneous presence of these may lead to cancellations. 
Eq.(ll) now gets replaced by 

(A-',„ + /4a^., + ^!3A-',3.) < 2.2 X 10-^' (^)' ^'-''^ 

With A'121 ~ 0.13, cancellation between the last two terms is unlikely as it requires A'131 ~ 2. 
The first two terms can cancel but the A'm is independently constrained from the neutrino 
mass. Its presence generates a large contribution to the electron neutrino mass induced 
through neutrino-gaugino mixing which can be given as, 

< i> >2 (6.27) 

msuSY 

The value of the induced sneutrino vev is sensitive to the MSSM parameters as we have 
seen above, but can be approximately written as. 



Requiring m^^ < 2eV leads for msusY ~ 100 GeV to 



Am < .04 
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It is seen that cancellations between the first two terms in eq.(12) are feasible and can 
allow A121 ~ 0.13 if this fine tuning is accepted. It must be added that the bound in the 
previous equation is quite sensitive to the MSSM parameters and for a large range in these 
parameters, the actual bound can be stronger than the generic bound displayed above. 

While wino and zino control the decay of the charm squark, its mass is mainly controlled 
by the large radiative corrections driven by the gluino mass. The unification of the gaugino 
mass parameters relates the two and leads to the above difficulty. Thus giving up this 
unification may open up a possibility of reconciling HERA events . Let us treat the gaugino 
masses Mi^2,3 at Mz as independent parameters. Then integration of the RG equation for 
the charm-squark from Mqut to Qq = 200 GeV leads to 



m 



2 (Qo) ~ ml (Mgut) + 0.77M| + 0.70M| + 0.024M2 



+ 1/2 M| COS 2/3 (1 -4/3 sin^6'iy) (6.29) 

If gaugino masses were to be unified at Mqut then M3 ~ 3.25M2 and Mi ~ O.5M2. Even 
in the absence of such unification, the physical gluino mass rrig ~ (1 + 4.2^)M3 must be 
greater than the charm squark mass if large \'i2i is to be avoided. This follows since in 
the converse case, the charm squark would predominantly decay to a gluino and a quark. 
This decay being governed by strong coupling, would dominate the other decays and would 
reduce B. The later is given in case of mg^ 3> rUg by 



32 vr V 0.1 

Such a tiny value of B would need unacceptably large \'i2i- It therefore follows that one 
must suppress the squark decay to gluino kinematically by requiring rUg > . Given this 
bound on M3 it follows from eq.(15) that 

M2<170GeV (6.30) 

if mg^ ~ 220 GeV. This bound on M2 is weaker than the one in the case of the gaugino mass 
unification, eq.(4). But it nevertheless cannot suppress the decay of squarks to chargino 
kinematically. It follows from Fig. 1 that one now approximately needs \'i2i > 0.08. This 
value is close to the 2a limit coming from the atomic parity violation but one would still 
need some cancellations to satisfy other constraints as discussed above. Thus giving up 
unification helps only partially. 

An alternative possibility is to allow for a -ve (mass)^ for the charm squark at the uni- 
fication scale. In view of the large positive contribution induced by the gluino mass such 
negative (mass)^ need not lead to colour breaking and may be consistent phenomenologi- 
cally. In fact a -ve (mass)^ for top squark has been considered in the literature in a different 
context. The universality is a simplifying feature of MSSM but it does not follow from any 
general principle. It does not hold in a large class of string based models which may allow 

®Note that fig.l is based on the assumption of Mi = O.5M2 but does not use any relation between M2 
and M3. 
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negative (mass) for some sfermions as well |25|. Such masses can also arise when SUSY is 



broken by an anomalous U{1) with some of the sparticles having -ve charge under this 
^(1). 

The large radiative corrections induced through the running in squark masses from a high 
~ Mgut to the weak scale has played an important role in this analysis. In contrast to the 
supergravity based models, this running is over a much smaller range in models with gauge 
mediated super symmetric breaking. But in these models, the initial value of the charm 
squark (mass)^ is positive and large with the result that these models are incompatible with 
the charm squark interpretation of HERA anomaly even without the radiative corrections 

The interpretation of HERA events in terms of stop may not suffer from the above 
mentioned difficulty encountered for the charm interpretation for two reasons. Firstly, the 
stop mass is reduced compared to the charm squark mass due to the possible large t^ — tji 
mixing as well due to the large top coupling. Secondly, this mass also involves one more 
parameter (the trilinear couplings A) compared to the charm squark mass. Thus while this 
is a less constrained possibility, imposition of the requirement that ~ 200 GeV would 
certainly lead to more a constrained parameter space than considered in model independent 
studies H^. 

In summary, we have shown that the charm squark interpretation of HERA events is 
possible only for large X'121 ~ 0(0. 1) in a large class of supersymmetric standard models 
characterized by a positive charm squark (mass)'^ at the GUT scale. The simplest and 
the most popular minimal supergravity model with universal boundary condition falls in 
this class. The required large value of R violating parameter is difficult to admit without 
postulating new physics and /or fine tuned cancellations due to constraints coming from the 
R-violating processes like neutrino mass, — > vr+z^P etc. However, the later data from 
HERA experiments did not substantiate their claim for a resonance. Treating the deviations 
at 200 GeV as statistical fluctuations, the HERA collaborations have now put constraints 
on the mass and generic couplings of leptoquarks. These bounds are now compatible with 
the results obtained from Tevatron |14|. 



^ See, K. Cheung, D. Dicus and B. Dutta |2T). 
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Conclusions 



The highly successful Standard Model has to be incorporated in larger theories to have a 
deeper understanding of Nature. Supersymmetric Standard Models have been built to solve 
the hierarchy problem which is encountered whenever Standard Model is being incorporated 
into a larger theory. Experimentally, one of the first signals for physics beyond Standard 
Model is provided by the atmospheric neutrino experiments which indicate the existence 
of atleast one non-zero neutrino mass. An interpretation of both solar and atmospheric 
neutrino anomalies interms of neutrino oscillations would lead to a specific pattern of the 
neutrino mass matrix. In this thesis, we studied the possibilities of generating such a pattern 
within the context of Supersymmetric Standard Models. 

Neutrino masses can be generated in the Supersymmetric Standard Model as it natu- 
rally allows for lepton number violation. Models which violate lepton number are classified 
into bilinear and trilinear lepton number violating models based on the dimensionality of 
the couplings in the superpotential. Firstly, we studied in detail models which violate lep- 
ton number only through bilinear (dimensionful) couplings. In these models, only three 
R-parity violating couplings determine the scale of the neutrino masses. The masses are 
naturally small due to Yukawa suppression if universal boundary conditions are assumed at 
the high scale. After presenting a general analytical structure of the neutrino mass matrix 
in these models, we have studied them in detail in the framework the Minimal Messen- 
ger Model (MMM) of gauge mediated supcrsymmctry breaking. The bilinear R-violating 
model within MMM, is very economical with essentially only four parameters, the three R- 
parity violating couplings and the supersymmetry breaking scale A, determining the three 
masses and three mixing angles of the neutrinos. We have shown that these models are 
too restrictive to provide simultaneous solutions for solar and atmospheric neutrino prob- 
lems. For simultaneous solutions with MSW conversion, we have numerically shown that 
these models do not provide the required mass hierarchy in the neutrino mass eigenvalues. 
On the other hand, for simultaneous solutions with vacuum oscillations, we have shown 
that these models cannot accommodate the two large mixings required. Thus, extensions 
of the minimal messenger model are inevitable for existence of simultaneous solutions for 
solar and atmospheric neutrino anomalies within these models. A simple extension which 
would allow for a negative /it parameter would be able to provide the required framework 
for simultaneous solutions. 

The other class of models we have considered are models with trilinear lepton number 
violation. These models, unlike the bilinear case, offer much larger freedom due to the 
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larger number of parameters present. One of the major results of this thesis is the effect of 
Renormalisation Group Evolution on the neutrino mass spectrum in the presence of only 
trilinear (dimensionless) lepton number violating couplings in the model. RG evolution 
generates non-zero soft bilinear lepton number violating couplings at the weak scale though 
these couplings are absent at the high scale. These would generate vacuum expectation val- 
ues (vevs) to the sneutrinos leading to non-zero mixing between ncutralinos and neutrinos. 
A 'tree level' neutrino mass is thus generated, which we call the 'RG induced tree level' 
mass. Analyses of the neutrino mass matrix in these models presented in the literature 
have been neglecting this important contribution which can drastically alter the neutrino 
mass spectrum in these models. We have analysed the neutrino mass spectrum taking into 
consideration RG induced tree level mass contribution in addition to the standard 1-loop 
induced masses within the framework of mSUGRA inspired MSSM. Our results show that 
the model favours vacuum solution to the solar neutrino problem while simultaneously al- 
lowing solutions for the atmospheric neutrino problem. However, there exist regions in the 
parameter space where the hierarchy in the neutrino mass eigenvalues is reduced due to a 
suppression of the tree level mass. In these regions, Large Angle MSW solution is preferred 
which is favoured by the latest data from super-Kamiokande. 

In addition to studying the neutrino mass spectrum in the R-parity violating scenarios, 
one can instead consider bounds on the lepton number violating couplings due to direct 
experimental limits on the neutrino masses. Since neutrinos attain majorana masses in 
this case, the most stringent limit on them comes from the non-observation of neutrinoless 
double beta decay. Using this experimental limit, bounds on the R-parity violating couplings 
can be derived. The existing limits do not consider the RG effects on the neutrino mass. 
But, as we mentioned above, the RG induced tree level mass would alter the neutrino mass 
spectrum drastically. We have included this effect and re-derived bounds on the lepton 
number violating couplings. In addition to this, we have also considered the effect of CKM 
matrix and discussed the basis dependence of these bounds. We have subsequently used 
these bounds in our discussion of the charm squark interpretation of the now defunct HERA 
anomalies. However, our analysis is of general nature and can be redone if a relevant case 
arises in the future colliders. 

General properties of neutrino masses in R-violating theories 

So far, in this thesis, we have concentrated on specific models of R-violation and super- 
symmetry breaking to study the structure of neutrino mass matrix. Though some of the 
properties of the neutrino mass matrix (essentially related to solutions of solar and atmo- 
spheric neutrino problems) depend on specific models of R-violation and type of supersym- 
metry breaking, most of the properties hold true in any of the R-violating scenarios. We 
summarise them as follows. 

(i) R -violating couplings 

Majorana masses for the neutrinos are produced in these models due to the presence of 
lepton number violating couplings. The neutrino masses are proportional to the square of 
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these couplings. In the hmit R-violating couplings are set to zero, there are no neutrino 
masses in these models and they represent the standard MSSM with R-parity. 

(ii) Hierarchical masses 

Neutrinos attain masses both at the tree level as well as at the 1-loop level. The tree 
level mass is generated indirectly through the bilinear soft L violating terms in the scalar 
potential. They are either present 'originally' in the scalar potential or generated at the 
weak scale through RG scaling. The total neutrino mass matrix is sum of the tree level 
and 1-loop contributions. At the tree level only one combination of the neutrinos attains 
mass. The other combinations attain mass at the 1-loop level only. The tree level mass is 
much larger compared to the 1-loop contribution for most of the R-violating theories. The 
hierarchy in the neutrino mass spectrum is typically characterised by . Since this 

factor is much smaller than one, large hierarchy in the neutrino mass spectrum is natural 
in these models. In some theories it is possible to find regions in parameter space where 
the tree level mass is suppressed and becomes comparable to the 1-loop level mass. 

(Hi) Large Mixing 

In most of these models, by suitable analytical approximations one can decouple the neutrino 
mass matrix into R-violation dependent part and R-violation independent part. The R- 
violation dependent part determines the mixing in most of the cases. The mixing angles 
are dependent on ratios of the R-violating parameters. When R-violating parameters are 
equal in magnitude large mixing is natural in the model Q. But this does not guarantee us 
that two large mixing angles can be simultaneously accommodated in these models. 

(iv) Role of Yukawas 

The R-violating parameters which determine the neutrino masses are always accompanied 
by down quark or charged lepton Yukawa couplings. Thus in the limit of vanishing charged 
lepton and down quark Yukawa couplings, neutrino masses also vanish in these theories (In 
models with only dimensionful R violation this is true only in the case of universal boundary 
conditions at the high scale). This dependence on the Yukawa couplings can be understood 
in terms of a U{1) symmetry. In the absence of down quark, charged lepton Yukawas and 
the n term, the superpotential is invariant under the U{1) with charges for the superfields: 

L = I ■,H2 = -I = -2 1; D"" = -I; U"" = I (7.1) 

and the rest of the superfields carry zero charge. I is an integer. The neutrino remains 
massless in this case. The Yukawas and the term break this symmetry and thus accompany 
the R-violating terms in the neutrino mass formulae. 

(v) Small R violation 



This may not hold true in some models with only A couplings. 
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Neutrinos are the only standard model fermions which attain mass through the soft 
sector in these theories. Whereas the loop induced masses are suppressed by the Yukawa 
couplings and are small, the tree level neutrino mass can be in general very large ~ 0{Mausy) 
in these models. These large neutrino masses have to be suppressed to be phenomenologi- 
cally meaningful. These masses can be suppressed either in a natural way or by a choice of 
suppression of the parameters: 

(a) By a suitable choice of the boundary conditions for bilinear lepton number violating soft 
terms at the high scale where supersymmetry is broken, one can bring down the neutrino 
masses to 0(100 MeV) in these models]^. For example, by choosing these terms to be zero 
at the high scale (as naturally happens in models with only dimensionless lepton number 
violation) or to be of the same magnitude as the other bilinear terms like at the high 
scale (as happens in models with universal boundary conditions) leads to bilinear lepton 
number violating soft terms which are only Yukawa suppressed at the weak scale due to 
RG scaling. This kind of suppression can generally lead to 0(100 MeV) neutrino masses. 

(b) One can always 'choose' the two contributions to the sneutrino vev to cancel each other. 
This naturally happens for some parameter space in mSUGRA inspired MSSM. But, in 
model independent analyses, which we discuss below the two parameters which contribute 
to the sneutrino vev or the sneutrino vev in general itself are free parameters of the model. 
In these cases, there is no natural choice of the sneutrino vev and it has to be chosen small 
in order to have a small neutrino mass. 

Both the above methods would not be able to suppress the neutrino mass sufficiently 
to give the right scale of ~ 0(eV). Thus one has to choose extremely small R- violation to 
generate 0(eV) neutrino masses in these models. This is possible as the neutrino masses are 
directly proportional to the R-violating couplings in these models. For example, to accom- 
modate solutions of solar and atmospheric neutrino problems, bilinear R-parity violating 
parameters are required to be 0{l{)~^)xMsusy Similarly the dimensionless A' couplings are 
typically required to be of 0(10^'^) in these models. Some models have been explored in 
literature where a natural way of such small R- violation can be achieved [Q, ^ 

(vi) Experimental Signatures 

R-violating theories allow for completely different experimental signatures in comparison 
to the standard MSSM. The additional L violating couplings have characteristic signatures 
some of which we have already seen in chapter 6. A comprehensive study of experimental 
signatures of R-violating theories has been recently presented in Q . The models presented 
in this thesis have restricted the R-violating parameter space as well as the standard super- 
symmetry soft parameters so as to provide simultaneous solutions to solar and atmospheric 
neutrino problems. Thus these models would have different mass spectrum and decay signa- 
tures. An analysis of this type combining neutrino masses and experimental signatures have 
been done in case of some models recently [^, ^. The other most important consequence of 
R-violation is that the LSP is no longer stable. One has to look for a new cold dark matter 

■^This may not generally hold true for some models with only trilinear A couplings. 



Conclusions 



102 



candidate in these models. 

Thus we see that R-parity violating theories provide a natural framework where small 
neutrino masses can be realised. In these models neutrino masses are calculable in terms 
of few basic parameters and thus can predict the existence of solutions for solar and atmo- 
spheric neutrino masses within specific models of supersymmetry breaking. Both bilinear 
and trilinear R-violating models have various interesting features associated with them 
making studies within these theories important and essential. 
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